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Abstract. We establish Schauder a priori estimates and regularity for solutions to a class of 
degenerate-elliptic linear second-order partial differential equations. Furthermore, given a C°°- 
smooth source function, we prove C°° -regularity of solutions up to the portion of the boundary 
where the operator is degenerate. Degenerate-elliptic operators of the kind described in our article 
appear in a diverse range of applications, including as generators of affine diffusion processes 
employed in stochastic volatility models in mathematical finance O [27] , generators of diffusion 
processes arising in mathematical biology O |TT] , and the study of porous media [S] [7] . 
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1. Introduction 

This article continues our development of regularity theory for solutions to the Dirichlet bound- 
ary value problem defined by a degenerate-elliptic operator. Degenerate-elliptic operators of the 
kind explored in our article can arise as generators of affine diffusion processes employed in sto- 
chastic volatility models in mathematical finance ^T7\ , generators of diffusion processes arising 
in mathematical biology [3l [11], and the analysis of porous media |6t [Tj, to name just a few 
applications. 

In [5] , in addition to other results, Daskalopoulos and Feehan obtained existence of solutions 
to a variational equation defined by the Heston operator [27] . We recall that the Heston operator 
serves as a useful paradigm for degenerate-elliptic operators arising in mathematical finance. In 
|17j . the authors proved global C"-regularity of solutions to the variational equation defined 
by the Heston operator, while in [19], the authors established Jif^ as well as C^'" and Cs'^^" 
regularity for those solutions, for all integers A; > 0. However, our Cs'" and Cs'^"*"" regularity 



results in |19j . although they provide an important stepping stone, are not optimal due to our 
reliance on variational methods. The purpose of the present article is prove analogues — for 
a broader class of degenerate-elliptic operators — of Schauder a priori estimates and regularity 
results for strictly elliptic operators in Chapter 6] . When coupled with results of [51 [TTl [T^] , 
we immediately obtain existence and Cs regularity for solutions to the Dirichlet boundary 
value problem, defined by a degenerate-elliptic operator, analogous to those expected from the 
Schauder approach for strictly elliptic operators in [26, Chapter 6]; uniqueness for broad class 
of linear second-order degenerate-elliptic operators, with the second-order (or Ventcel) boundary 
conditions of the kind implied by our choice of Daskalopoulos-Hamilton Cg Holder spaces |6], 
is a consequence of the weak maximum principle discussed by the first author in |16]. 

To describe our results in more detail, suppose ^ ^ H is a domain (possibly unbounded) in 
the open upper half-space HI := W'-~^ x M+, where d > 2 and M+ := (0, oo), and di^ := Off nM 
is the portion of the boundary 80' of which lies in H, and Oq^ is the interior of dM D 80, 
where 811 = W^'^ x {0} is the boundary of M := M'^"^ x M+ and M+ := [0, cx)). We assume 8o0 
is non-empty and consider a linear second-order elliptic differential operator, L, on which is 



degenerate along 9o^. In this article, when the operator L is given by (1.3), we prove an a priori 
interior Schauder estimate and higher-order Holder regularity up to the boun dary portion, 80^ 
— as measured by certain weighted Holder spaces, Cs'^'''"(^) (Definition 2.3) — for solutions to 
the elliptic boundary value problem, 

Lu = f on0, (1.1) 
u = g on8i0, (1.2) 

where / : 1^ — )• M is a source function and the function g : 8i0' — )• M prescribes a Dirichlet 
boundary condition. We denote := 0'U8o0' throughout our article, while = G\J8G denotes 
the usual topological closure of & in . Furthermore, when / G C°^(^), we will also show that 
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u £ C°°{^) (see Corollary 1.9). Since L becomes degenerate along 9o^) such regularity results 
do not follow from the standard theory for strictly elliptic differential operators [26. i30j. 

Because the coefficient, b'^, will be assumed to obey a positive lower bound along doff, no 
boundary condition is prescribed for the equation (1.1) along doff. Indeed, one expects from [6] 
that the problem (1.1), (1.2) should be well-posed, given / G Cg{ff) and g G C(9i^) obeying 

). The degenerate- 



mild pointwise growth conditions, when we seek solutions in C^^"(^) H Cio 
elliptic operator considered in this article has the forrr|^ 

Lv:=-Xdtr{aD'^v)-b-Dv + cv on ^, u G C°°(^), (1.3) 

where x = (xi, . . . , Xd) are the standard coordinates on W^. Occasionally we shall also need 

Lov:= {L - c)v = -XdtT{aD^v) - b ■ Dv on ff, veC^{ff). (1.4) 



Throughout this article, we assume that the coefficient functions a, 6, c of L (and Lq) are defined 
on ff, the matrix (a*-') is symmetricj^ and there is a positive constant bo such that 

b"^ > bo on doff. 



In 



We shall call L in ( 1.3 ) an operator with constant coefficients if the coefficients a, b, c are constant. 
[T9] . we proved existence and uniqueness of a solution, u G Cg~^'^°{ff) (1 C{ff) for some 
ao G (0,1), to (1.1), (1.2) when diff obeys a uniform exterior cone condition with cone 
K, and L is the ehiptic Heston operator, and / G C°°(^) n Cb{ff) and g G C°°{Aff). (The 



Holder exponent, oq, depends on the coefficients of L and the cone K.) In ^1.1, we state the 
main results of our article and set them in context in ^1.2 where we provide a survey of previous 
related research by other authors. In §1.3[ we indicate some extensions of methods and results 



in our article which we plan to develop in subsequent articles. We provide a guide in { 1.4 to the 
remainder of this article and point out some of the mathematical difficulties and issues of broader 
interest. We refer the reader to 31.51 for our notational conventions. 



1.1. Summary of main results. We summarize our main results. Here, our use of the term 
"interior" is in the sense intended by p], for example, U C ^ is an interior subdomain of a 
domain ff if U G ff and by "interior regularity" of a function u on ff, we mean regularity 
of n up to 5o ^ — see Figure |1.1[ 

Our first main result is the following analogue of |Sl Theorem 1.1.3] (for a related degenerate- 
parabolic operator ( |1.25 ) and d = 2), ^ Theorem 3.1] (for a related degenerate-parabolic operator 
(1.25) with d > 2), and [26, Corollary 6.3 and Problem 6.1] (strictly elliptic operator). We refer 



the reader to Definitions 2.1, 2.2 and 2.3 for descriptions of the Daskalopoulos-Hamilton family 
of Cs and Cg Holder norms and Banach spaces. For any [/ ^ H, we denote 



-,k,a,fT\ and ll&ll^fc,a^r"r\ := ^ ^ 1 1 6' ||^fc,t:t 



1=1 



(1.5) 



Theorem 1.1 (A priori interior Schauder estimate). For any a G (0,1), integer k > 0, and 
positive constants bo, Ao, do, A, u, there is a positive constant, C = C{bo,d,do,k,a,Xo,A.,i'), 



^The operator —L is the generator of a degenerate-diffusion process with killing. 

9 

The assumption of symmetry is just for convenience when applying changes of variables and is easily obtained 



by replacing a'-' with a*- 
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Figure 1.1. Boundaries and regions in Theorem 1 1 . 1 1 and Remark 1.2 



such that the following holds. Suppose^height^ff) < v and the coefficients a,b,c of L in (1.3) 



belong to Cs'°^{ff) and obey 



K,0>Ao|^|' on 

b'^ > bo on di 



If u G Cs''^~^°'{^) and 0"' G ^ is a subdomain such that d.\si{di0'' ,di&) > do, then 



\u\\^k,2 + a:^^,~^ < C 



u 



C{6) 



(1.6) 
(1.7) 
(1.8) 

(1.9) 



Remark 1.2 (A priori interior Schauder estimate). The case where k = and the domain is a 
half-bah, ^ = i3+(x''), the coefficients, a,b, of L are constant, c = 0, and u G C°°(5+(x'')) is 
given by Corohary 6.8 Theorem 3.2 relaxes those conditions to allow u G Cg~^°' {Bj^^{x^)) and 
arbitrary c G M; Theorem 8.1 fur tner relaxes the conditions on L to allow for variable coefficients, 

relaxes the constraint A; = to allow for arbitrary integers 



a,b,c, in Cf{Bj {x )); Theorem 



1.1 



A; > 0; finally. Theorem 
allow for arbitrary domains 



8.3 



is proved in ^8.1 



where we relax the constraint that = (x^) and 



It is considerably more difficult to prove a global a priori estimate for a solution, n G 



k,2+a 



when the intersection do^Cidi^ is non-empty and we do not consider that problem in this article, 
but refer the reader to [HI §1.3] for a discussion of this issue. However, the global estimate in 
Corollary 1.3 has useful applications when di^ does not meet do( 
define the strip, 

S := 



For a constant > 0, we 



' X (0,1.), 



(1.10) 



■^If we had allowed heiglit(i^7') = oo, we would need to modify our definition of Holder norms to provide a weight 
for additional control when Xd ^ oo because the coefficient matrix, Xda, for D^u would be unbounded due to 



( 1.7l. Weighted Holder norms of this type were used by the authors in [18], for this reason, for the corresponding 
parabolic operator, —dt + L, on (0, T) x H. 
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and note that OqS 



t>d-l 



X {0} and diS 



Corollary 1.3 (A priori global Schauder estimate on a strip). For any a G (0, 1), positive 
constants Xq, bo, A, v, and integer k>0, there is a positive constant, C = C{k,a,i',d, Xo,bo, A), 



S 



such that the following holds. Suppose the coefficients of L in (1.3) belong to Cs'"^{S), where 
'^^^ X (Ojiy) as in (1.10), and obey 

(1.11) 
(1.12) 
(1.13) 



If u G Cs''^~^°'{S) and u 



and, when c>0 on S, 



< A, 

«,0>Ao|e|' onS, V^GM^ 
b'^ > bo on doS. 

on diS, then 



u 



< C 



+ u 



C(5) 



< C\\Lu\ 



(1.14) 



(1.15) 



Remark 1.4 (A priori global Schauder estimate on a strip). For an operator, L, with constant 



coefficients, o, b, c, an a priori global Schauder estimate on a strip is proved as Corollary 7.2 



The Green's function for a operator L in ( |1.3[ ) with constant coefficients can be extracted from 
Appendix [B| where we construct explicit solutions to Lu = / on HI and prove the following 
elliptic analogue of the existence result [6l Theorem 1.1.2] for the initial value problem for a 
degenerate-parabolic model (1.25) on a half-space for the linearization of the porous medium 



equation (1.24) 



Theorem 1.5 (Existence and uniqueness of a C°°{M) solution on the half-space when L has 
constant coefficients). Let L be an operator of the form (1.3) and require that the coefficients, 
a,b,c, are constant with b'^ > and c > 0. 
u£C^{U), to Lu = f onM. 



iff^Q 



then there is a unique solution, 



Again, it is considerably more difficult to prove existence of a solution, u, in Cs' ' "{(7) or 
n C{&), to (1.1), (1.2) when the intersection do^ H is non-empty. We do not 
consider that problem in this article either and again refer the reader to [19, §1.3] for a discussion 
of this issue. However, in the case of a strip, dii? does not meet do^ and we have an existence 
result. Theorem 1.6, for an operator with variable coefficients. In §1.3[ we discuss additional 
existence results which should also follow from Theorems LI and L5 when d^i? is curved and 
did 



is empty. 



tk,2+a 



(S) solution on a strip S). Let a G (0, 1) 



Theorem 1.6 (Existence and uniqueness of a C, 

^ci-i ^ (0, zy) be as in (1.10), and let k > be an integer. Let L be an operator 

tk,a 



let u > and S 



as in (1.3). If f and the coefficients of L in (1.3) belong to C, 



for some positive constants bo, Xo, then there is a unique solution, u ^ C. 
value problem, 

Lu = f on S, 
u = on diS. 



(S) and obey ( |LT2| ) and (L13 ) 

A;,2+Q!^^-j^ to the boundary 



(1.16) 
(1.17) 



Remark 1.7 (Existence and uniqueness of a solution on a strip). For an operator, L, with constant 



coefficients, o, 6, c, existence and uniqueness of a solution on a strip is proved as Corollary B.4 
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The preceding existence and uniqueness result on a strip leads to the following analogue of |26| 
Theorem 6.17] and is proved in j ]8.2[ 

Theorem 1.8 (Interior Cs '^"'""-regularity). For any a G (0,1) and integer k > 0, the following 
holds. Assume that the coefficients of L in (1.3) belong to Cs'"'{^) and obey (1.7) and (1.8) for 
some positive constants bo,Xo. If u G C^(^) o6e?/s|^ 



u G Ci 



x^D u G 



XdD^u = 



, and 
on doiff, 



Lu G C: 



(1.18) 
(1.19) 



then u ^ C 



k,2+a 



Given Theorem |1.8[ one immediately obtains the following degenerate-elliptic analogue of 
the C°°-reguIarity result for the degenerate-parabolic model for the linearization of the porous 
medium equation [6l Theorem 1. 1.1]. 



Corollary 1.9 (Interior C°°-regularity). Assume that the coefficients of L in (1.3) belong to 
C°°(^) and obey (1.7) and (1.8) for some positive constants 5o,Ao. If u ^ C'^{iff) obeys (1.18) 



for every integer k >0, so Lu G C°°{^), together with (1.19), then u G 

Remark 1.10 (Regularity up to the "non-degenerate boundary"). Regarding the conclusion of 
Theorem |1.8[ standard elliptic regularity results for linear, second-order, strictly elliptic operators 
[26} Theorems 6.19] also imply, when k > 0, that u G C''^^'' 
with / G C'='"( 



U diff) if u solves pTTl), (fO) 



'U Oi^) and g G C'=+2'°(^ U 9i^), and di^ is C7'=+2.". Because our focus in this 
article is on regularity of u up to the "degenerate boundary", d^ff, we shall omit further mention 
of such straightforward generalizations. 



Finally, we refine our existence results in [T^ when d = 2 for the Heston operator, 
Av := ivm.x, + 2gavm,x', + cr^'v^^^r.^) - (cq - q - X2/2)vxi - k{0 - X2)vx2 + cqv, (1.20) 



^2 / „ 2 



where g > 0, cq > 0, k > 0, > 0, o" > 0, and g G ( — 1, 1) are constants (their financial interpreta- 
tion is provided in [27|), and v G C°^{M). In particular, we give analogues of the existence results 
|26^ Theorems 6.13 & 6.19] for the case of the Dirichlet boundary value problem for a strictly 
elliptic operator. 



Theorem 1.11 (Existence and uniqueness of a Cs''^^" solution to a Dirichlet boundary value 
problem for the Heston operator). Let a G (0, 1) and let k > be an integer, let K be a finite 
right- circular cone, and require that d\G obeys a uniform exterior cone condition with cone K. 



If f & Cs'"{0nCbiff) and 



Co > if heig: 
Co > if heig: 



= oo, 
< oo. 



(1.21) 



then there is a unique solution, 

u G c^2+"(^) n C{& u ai^) n Cfe(^), 

to the boundary value problem for the Heston operator, 

Au = f on 
u = ondi&. 



(1.22) 
(1.23) 



We write Du, XdD^u £ C{ff) as an abbreviation for u^^ , XdUx^xj G C{tf), for 1 < i, j < d and write XdD^u = 



on doff as an abbreviation for limg 



3,0 XdD^u{x) = for all 2:" G doff- 
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Remark 1.12 (Schauder a priori estimates and approach to existence of solutions). As we explain 
in [m §1.3], the proof of existence of solutions, u G Cs''^~^°'{^) nC{S'), to the boundary value 
proble m, ([lT| ), ([OJ), given / G Cs'"{0 and g G C(^), is considerably more difficult when 
do^' n di0' is non-empty because, unlike in [6], one must consider a priori Schauder estimates 
and regularity near the "corner" points of the subdomain, ^ C M, where the "non-degenerate 
boundary", di0', meets the "degenerate boundary", do0'. 



Given an additional geometric hypothesis on ^ near points in Oq^ H d\G ^ the property that 
u G C(i^U d\G^ n Ch{^) in the conclusion of Theorem 1 1 . 1 1 1 simplifies to u G 



Corollary 1.13 (Existence and uniqueness of a globally continuous Cs solution to a Dirichlet 
boun dary value problem for the Heston operator). // in addition to the hypotheses of Theorem 
\l.ll\ the domain, ff, satisfies a uniform exterior and interior cone condition on Oq^? D di^ with 
cone K in the sense of [19j, then u G Cs'^^"(^) H 



Remark 1.14 (Existence of solutions to a Dirichlet boundary value problem). Theorem 1.11 and 
Corollary 1.13 should generalize to from the Heston operator ^ in (1.20) to an operator L in (1.3) 



with coefficients and d> 2. 

1.2. Survey of previous related research. We provide a brief survey of some related research 
by other authors on Schauder a priori estimates and regularity theory for solutions to degenerate- 
elliptic and degenerate-parabolic partial differential equations most closely related to the results 
described in our article. 



The principal feature which distinguishes the equation (1.1), when the operator L is given by 



(1.3), from the linear, second-order, strictly elliptic operators in and their boundary value 
problems, is the degeneracy of L due to the factor, Xd, in the coefficient matrix for D^u and. 



because 6o > in (1.3), the fact that boundary conditions may be omitted along x^ = when 
we seek solutions, n, with sufficient regularity up to Xd = 0. 

The literature on degenerate elliptic and parabolic equations is vast, with the well-known 
articles of Fabes, Kenig, and Serapioni jl21ll3|. Fichera [HI [22], Kohn and Nirenberg [29], Murthy 
and Stampacchia |33| l34j and the monographs of Levendorskii [31j and Olemik and Radkevic 
[351 IMl EZj, being merely the tip of the iceberg. 

As far as the authors can tell, however, there has been relatively little prior work on a priori 
Schauder estimates and higher-order Holder regularity of solutions up to the portion of the domain 
boundary where the operator becomes degenerate. In this context, the work of Daskalopoulos, 
Hamilton, and Rhee [6l[71|38] and of Koch stands out in recent years because of their introduction 
of the cycloidal metric on the upper-half space, weighted Holder norms, and weighted Sobolev 
norms which provide the key ingredients required to unlock the existence, uniqueness, and higher- 



order regularity theory for solutions to the porous medium equation (1.24) and the degenerate- 



parabolic model equation (1.25) on the upper half-space given by the linearization of the porous 
medium equation in suitable coordinates. 

Daskalopoulos and Hamilton ^ proved existence and uniqueness of C°° solutions, u, to the 
Cauchy problem for the porous medium equation [6l p. 899] (when d = 2), 

d 

-ut + Y,iu")-.x.=0 on(0,r)xM^ u{-,0)=g on (1.24) 

i=l 

with constant m > 1 and initial data, g > 0, compactly supported in R*^, together with C°°- 
regularity of its free boundary, d{u > 0}, provided the initial pressure function is non-degenerate 
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(that is, Du^~^ > a > 0) on boundary of its support at t = 0. Their analysis is based on their 
development of existence, uniqueness, and regularity results for the linearization of the porous 
medium equation near the free boundary and, in particular, their model linear degenerate operator 
[U p. 901] (generalized from d = 2 in their article). 



Lu 



^Ux„ u G C°°(]H), 



(1.25) 



where /3 is a positive constant, analogous to the combination of parameters, 2k0/cj^, in (1.20), 
following a suitable change of coordinates [H p. 941]. 

The same model linear degenerate operator (for (i > 2), was studied independently by Koch 
[28} Equation (4.43)] and, in a Habilitation thesis, he obtained existence, uniqueness, and regu- 
larity results for solutions to (1.24) which complement those of Daskalopoulos and Hamilton [6]. 



Koch employs weighted Sobolev space methods, Moser iteration, and pointwise estimates for the 
fundamental solution. However, by adapting the approach of Daskalopoulos and Hamilton [6], 
we avoid having to rely on difficult pointwise estimates for the Green's function for the operator 
L in (1.3). Although tantalizingly explicit — see P [HI HSl [201 [27] for the Green's function and 



fundamental solution of the elliptic and parabolic Heston operator (1.20) and Appendix B 



these kernel functions appear quite intractable for the analysis required to emulate the role of 
potential theory for the Laplace operator in the traditional development of Schauder theory in 

While the Daskalopoulos-Hamilton Schauder theory for degenerate-parabolic operators has 
been adopted so far by relatively few other researchers, it has also been employed by De Simone, 
Giacomelli, Kniipfer, and Otto in |8t 125^ IM] and by Epstein and Mazzeo in jlOj . 



1.3. Extensions and future work. We defer to a subsequent article the development of a 

estimates, existence, and regularity theory for solutions u to 



priori global Schauder C. 
the elliptic boundary value problem (1.1), (1.2) when / and the coefficients, a,b,c, of L in (1.3) 
belong to Cs'"{^), the boundary data function g belongs to Cs''^~^'^{^), and has boundary 
portion diff of class (7^^+^'" and C'^'^"'"'^-transverse to do^. For reasons we summarize in |19[ §1.3], 
the development of global Schauder a priori estimates, regularity, and existence theory appears 
very difficult when the intersection Oq^ D di^ is non-empty. 

However, if ^ C M'^ is a bounded domain and L is an elliptic linear second-order partial 



1.1 



differential operator which is equivalent to an o perator L^" of the form (1.3) in local coordinates 

will quickly lead to a global Cs''^^°' a priori 
Moreover, for g € Cs'"^^"^' 



near every point x G (9(^, then Theorem 

estimate for u if d& is of class Cs''^^°. Moreover, for g € Cs'^~^"{^), the method of proof of 
|30i Theorem 6.5.3] (or indeed ^ Theorem II. 1.1]) should adapt to give existence of a solution 



u G c^'^+'^i^) to (ITT]), ([m 



As we noted in Remark 1.14, we expect our existence results (Theorem 1.11 and Corollary 



1.11) for solutions to the Dirichlet boundary value problem (1.22), (1.23) to generalize from the 



case of the Heston operator. A, in (1.20) on subdomains of the half-plane to operators of the form 
L in (1.3) on subdomains of the half-space with Cs'"'^^) coefficients, a, b, c. These generalizations 



may be developed in two ways. First, the proof of Theorem 1.11 relies on existence and regularity 



theory for solutions to a variational equation defined by (1.22) and a choice of suitable weighted 



Sobolev spaces in [21 [13 [TH] ; we expect that analysis to extend without difficulty to operators of 



the form L in (1.3). Second, we expect the a priori interior Schauder estimates that we develop 



in this article, which are in the style of [26t Corollary 6.3], to extend to more refined and sharper 
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a priori "global" interior Schauder estimates, in the style of [26, Theorem 6.2, Lemmas 6.20 
and 6.21]. Aside from facilitating "rearrangement arguments", we expect such "global" a priori 
interior Schauder estimates — relying on a choice of suitable weighted Holder spaces similar to 
those employed in |26t Chapter 6] — to permit the use of the continuity method to prove existence 
of solutions within a self-contained Schauder framework and this theme will be developed by the 
authors in a subsequent article. 

While our a priori Schauder estimates rely on the specific form of the degeneracy factor, Xd, of 



the operator L in (1.3) on a subdomain of the half-space, we obtained weak and strong maximum 
principles for a much broader class of degenerate operators in [16j. Therefore, we plan to extend 
the a priori Schauder estimates and regularity theory for degenerate-elliptic operators such as 

Lv = -'di^:[aD'^v)-b-Dv + cv on ^, ?; G C°°(^), 

where "d G C^^{^) and -i? > on a subdomain 0" <Z with non-empty boundary portion 
do& = int({x G d& : = 0}). 



1.4. Outline and mathematical highlights of the article. For the convenience of the reader, 
we provide a brief outline of the article. In ^ , we review the construction of the Daskalopoulos- 
Hamilton- Holder families of norms and Banach spaces [6j. 

In ^ we derive a priori local estimates for derivatives of solutions, u, to Lu = on half- 
balls B^^{xo) C IH centered at points G dM, when L has constant coefficients. However, our 
method of proof differs significantly from that of Daskalopoulos and Hamilton [6|, who apply 
a comparison principle for a certain non-linear parabolic operator and which directly uses the 
fact that this operator is parabolic. We were not able to replace their "parabolic" comparison 
argument by one which is suitable for the elliptic operators we consider in this article. Instead, 
we employ a simpler approach using a version of Brandt's finite-difference method [1] to estimate 
derivatives in directions parallel to dM and methods of ordinary differential equations to estimate 
derivatives in directions normal to dM. 

In Q we adapt and slightly streamline the arguments of Daskalopoulos and Hamilton in [6] for 
their model degenerate-parabolic operator (1.25) to the case of our degenerate-elliptic operator 



(1.3) and derive a C'' a priori estimate of the remainder of the first-order Taylor polynomial of a 

B+{xo). 



function n on a half-ball 

In ^ we obtain a priori local interior Schauder estimates for a function ii on a ball B, 



G 



and dM, again when L has 



where we keep track of the distance between the ball center, x^ 
constant coefficients. 

In ^ we apply the results of the previous sections to prove our main Cg~^°' a priori interior 
local Schauder estimate (Theorem 3.2) for an operator L with constant coefficients on a half-ball. 

In ^ we prove Cs '^^" a priori interior local Schauder estimate (Theorem 7.1) and a global a 
priori global Schauder estimate on a strip (Corollary 7.2), both when L has constant coefficients. 



Injsj we relax the assumption in the preceding sections that the coefficients of the op erato r L 
in (1.3) are constant and prove a C^^" a priori interior local Schauder estimate (Theorem 8.1) for 
a function, n, on a half-ball, (x'^) when L has variable coefficients. We then prove a Cs' "''"a 



of Theorem 



Theorem 



the G 



1.1 



1.6 



and Corollary 



1.3 



priori local inter ior Schauder estim ate for arbitrary /c G N (Theorem 8.3) and complete the proofs 

Next, we prove our global Cs''^~^°'{S) existence result on st rips, 
and complete the proofs of our m ain Cs'^"*"" re gulari ty result. Theorem 



1.8 



existence results. Theorem 



1.11 



and Corollary 



1.13 



and 



for solutions to a Dirichlet 



boundary value problem for the Heston operator. 
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We collect some additional useful results and their proofs in several appendices to this article. 
In Appendix [Aj we prove a comparison principle for operators which include those of the form 
L in (1.3) with c > (rather than c > cq for a positive constant cq) when the domain, ^, is 
unbounded but has finite height, extending one of the comparison principles in |16j . In Appendix 
[B| we prove Theorem 1.5 In Appendix [Cj we summarize the interpolation inequalities and 
boundary properties of functions in weighted Holder spaces proved in [6j and [18] . 



1.5. Notation and conventions. In the definition and naming of function spaces, including 
spaces of continuous functions and Holder spaces, we follow Adams [2j and alert the reader to 
occasional differences in definitions between [2] and standard references such as Gilbarg and 
Trudinger [26j or Krylov |30j . 

We let N := {0, 1, 2,3,.. .} denote the set of non-negative integers. If S" C M"', we let S denote 
its closure with respect to the Euclidean topology and denote dS := S\S. For r > and x*^ S W^, 
we let Br{x^) := {x G M*^ : |x — x"! < r} denote the open ball with center and radius r. We 
denote B+{x^) := Br{x^) Pi BI when x° G 91HI. When x° is the origin, O £ R"^, we denote Br{x^) 
and by Bj- and B^ for brevity. 

If V C U C M"' are open subsets, we write V U when U is bounded with closure U C V. By 
suppC, for any ( G C(M'^), we mean the closure in M'^ of the set of points where C 7^ 0- 

We use C = C{*, ...,*) to denote a constant which depends at most on the quantities appearing 
on the parentheses. In a given context, a constant denoted by C may have different values 
depending on the same set of arguments and may increase from one inequality to the next. 



2. Preliminaries 

In this section, we review the construction of the Daskalopoulos-Hamilton-Holder families of 
norms and Banach spaces [6j. 

We first recall the definition of the cydoidal distance function, s (•,•), on H by 

I 1 2 1 

s{x\x^) := \^ -X \ y^l^^2 g jj^ (2.1) 

\/^d + + — x'^l 

where x* = {x\, . . . , x^^), for z = 1,2, and \x^ — denotes the usual Euclidean distance be- 
tween points G M'^. Analogues of the cycloidal distance function (2.1) between points 
{t^,x^), (i^,x^) G [0,oo) X H, in the context of parabolic differential equations, were introduced 
by Daskalopoulos and Hamilton in [6l p. 901] and Koch in [28, p. 11] for the study of the porous 
medium equation. 

Observe that, by (2.1), 

s(x,x°) < |x-xY^^ Vx,x°G]HI. (2.2) 

The reverse inequality take its simplest form when x*^ G dM, so x^ = 0, in which case the 
inequalities Xd < |x — x'^l and 

|x — x^'l = s(x, x^)^/xd^+~\x~^^x^\ < s{x, x'^) \/ 2\x — x^\, 

give 

|x-x°| < 2s(x,x°)^ Vx G M, x° G 5M. (2.3) 

Following [21 §1.26], for a domain U C H, we let C{U) denote the vector space of continuous 
functions on U and let C{U) denote the Banach space of functions in C{U) which are bounded 



REGULARITY OF SOLUTIONS TO DEGENERATE-ELLIPTIC EQUATIONS 11 

and uniformly continuous on [/, and thus have unique bounded, continuous extensions to [/, with 
norm 

Il^llc(t7) ■= sup|u|. 
U 

Noting that U may be unbounded, we let Cioc(f^) denote the linear subspace of functions u £ C{U) 
such that u € C{V) for every precompact open subset V U. We let Ci,{U) := C{U) fl L°°{U). 
Daskalopoulos and Hamilton provide the 

Definition 2.1 (C" norm and Banach space). |6j p. 901] Given a G (0, 1) and a domain [/ C H, 
we say that u G Cg{U) if u G C{U) and 

where 

IMc^iu) ■= Mcrm + \Mc{u)^ (2-4) 

and 

\u{x^) — u{x'^)\ 



^""^^'^^^^ JXu six\x^r ■ ^'-'^ 

We say that u G C°{U) if u G Cg{V) for all precompact open subsets V U^, recalling that 
U := U U doU. We let Cf^^^{U) denote the linear subspace of functions u G Cg{U) such that 
u G Cg{V) for every precompact open subset V <^U . 



It is known that Cf{U) is a Banach space [HI §1.1] with respect to the norm (2.4). 

We shall need the following higher-order weighted Holder and Cs norms and Banach 
spaces pioneered by Daskalopoulos and Hamilton [6]. We record their definition here for later 
reference. 

Definition 2.2 (C, " norms and Banach spaces). [6l p. 902] Given an integer A: > 0, a G (0, 1), 
and a domain [/ C H, we say that u G Cs '"(f/) if u G C^{U) and 



where 



E WD^^Wcm^ (2-6) 

where ^ := (A, . . . , /3rf) G N'^ and 

ll'"llcf(C7) ■= ll"llc(l/) + Mcf{C7)- 
When A; = 0, we denote C°'"(C7) = Cf{tJ). 

Definition 2.3 (Cr^+" norms and Banach spaces). [U pp. 901-902] Given an integer A: > 0, 
a constant a G (0,1), and a domain [/ C H, we say that u G Cs'^'''"(C7) if u G Cs^"'^'"(L^), 
the derivatives, -D^u, /3 G N'^ with = /c + 2, of order A; + 2 are continuous on U , and the 
functions, XdD^u, /3 G N'^ with |/3| = A; + 2, extend continuously up to the boundary, dU, and 
those extensions belong to Cf{JJ). We define 

\\u\\^k,2+c^jj^ := ||u||^fc+i,a^^^ + ^ WxdD'^uWc^^^ijy (2.7) 

\0\=k+2 
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We say thalj^n G Cs'^'^°(C/) if u G Cs''^~^°'{V) for all precompact open subsets V (£ U_. When 
fc = 0, we denote C°'^+"(t/) = C2+"(C7). 

For any non- negative integer k, we let Cq{U) denote the linear subspace of functions u G C^{U) 
such that u G C^{V) for every precompact open subset V <^U_ and define C^{U) := r]k>oCQ{U_). 
Note that we also have C^iU) = nfc>oCs'°(I7) = nk>oCs'^'^"{U). 

3. Interior local estimates of derivatives 

As in [6j, we begin with the derivation of local estimates of derivatives of solutions on half- 
balls B^{xq) centered at points G dM, but the method of proof differs significantly from 
the method of proof in [U §1.4 & 1.5]. In [6], Daskalopoulos and Hamilton apply a comparison 
principle to a suitably chosen function, defined in terms of the derivatives (see the definitions of 
Y at the beginning of [6i §1.5] and of X in the proof of [6j, Corollary 1.5.3]). Their comparison 
principle directly uses the fact that the operator is parabolic, and we were not able to replace the 
"parabolic" comparison argument by one which is suitable for the elliptic operators we consider in 
this article. (The Daskalopoulos-Hamilton approach can be viewed as a variant of the Bernstein 
method — see the proof [30, Theorem 8.4.4] in the case of the heat operator and [301 Theorem 
2.5.2] in the case of the Laplace operator.) 

Instead, we apply a combination of finite-difference arguments, methods of ordinary differential 



equations, and, in this section, restrict to the homogeneous of version of the equation (1.1) with 
/ = 0. We adapt Brandt's finite-difference method [1] (see also [261 §3.4]) to obtain a priori local 
estimates for D^u, where /3 G N*^ is any multi-index with non-negative integer entries of the form 
/3 = (/?!,..., /3'^-\ 0). The method of Brandt also uses a comparison principle, but it is applied 
to finite differences, instead of functions of derivatives of u, such as X and y in ^ §1-5]. Brandt's 
approach is also mentioned by Gilbarg and Trudinger in [261 p. 47] as an alternative to the usual 
methods for proving a priori interior Schauder estimates such as [26^, Corollary 6.3]. We are able 
to apply the finite-difference estimates method not only on balls Br{x^) d EI as in 0], but also on 
half-balls i?+(xo) C EI centered at points x^ G dM because the degeneracy of the elliptic operator 



L in (1.3) along doB:^{xQ) and the fact that b'^ > along doB^{xo) (see (3.2)) implies that no 
boundary condition need be imposed along doB:^{xo). 

In §3.1| we summarize the inte rior local Schauder estimate and regularity results we will prove 



in sections pi ^ androl In ^3.2, we develop interior local estimates for derivatives D^u when 



f^d = and in §3.3[ we extend those estimates to case /3rf > 0. 



3.1. A priori interior local Schauder estimate and regularity statements in the case 
of constant coefficients. Throughout sections^^^ and\^ we further assume the 

Hypothesis 3.1 (Constant coefficients and positivity). The coefficients, a,b,c, of the operator 
L in (1.3) are constant; there is a positive constant, Aq, such that|^ 

«,0>Ao|eP, V,^GM'^; (3.1) 

andlH 

b'^ = bo> 0. (3.2) 



^In [6l pp. 901-902], when defining the spaces (^^'"(jz/) and Cs''^^"{ji/), it is assumed that .(2/ is a compact 
subset of the closed upper half-space, H. 



Condition (|3.1l) is first used in the proof of Lemma 



5.1 



Condition (|3.2[) is required by our maximum principle (Lemma|A.lJand CoroIlary|A.2|. Our maximum principle 



in turn required in ^ sections jZ] and [s] depend on and sections jsj [t] and js] each depend on sections sjandjsj 
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The condition (3.2) is first required in the proof of Lemma 5.1 When the coefficients of L are 



constant, we denote 



A 



+ max 16*1 + IcI. 



(3.3) 



Our main goal in sections [3j |4j [5} and [6] is to prove the following versions of Theorems 1.1 and 1.8 



when A; = and L has constant coefficients and the domain is a half-ball, B^^^{x^) with G (9H. 

Theorem 3.2 (A priori interior local Schauder estimate when L has constant coefficients). For 

any a G (0, 1) and constants r and tq with < r < tq, there is a positive constant, C = 
C(a, r, ro, d, Ao, 6o, A), such that the following holds. If ^ dM and u Cg~^°'{Bj^^^{x^)), then 



< C 



(3.4) 



Our goal in the remainder of this section is to derive a priori estimates for Du and XdD^u on 



half-balls, B^. 



centered at points x G dM. Because our operator, L, is invariant with respect 



to translations in the variables (xi, . . . ,Xd-i) when the coefficients, a,b,c, are constant, we can 

3+ 

ro 



assume without loss of generality that x*^ is the origin, O G M*^, and write B^ {x'^^ 



B^{x^) = B^ in our proof of Theorem 



3.2 



i3+ and 



3.2. Interior local estimates for derivatives parallel to the degenerate boundary. To 

derive a priori local estimates for D^^u, for /? G N'^ with j3d = 0, it will be useful to consider the 
following transformation, 

u{x) =: v{y), X G H, 

where y = (f){x) := x + ^Xd and = {^i, ■ ■ ■ , ^d-i,S,d) S M*^. We choose such that 

:=-676^ Wi^d, (d = 0, 

where we have used assumption (3.2) that b'^ > 0. Note that (/> is a diffeomorphism on 
restricts to the identity map on dM. We now consider the operator Lq defined by 

Lqu{x) =: Lov{y), x G H, 

and by direct calculations we obtain 



(3.5) 

(3.6) 

which 



Lqv = -yddJ-^Vy^y^ - b'vy 



on 



(3.7) 



where 



-id 



a''^ + ^io'^'^, yi / d, 

dd 



b' + vi ^ d. 



+ iii,a''\ yij^d, 



(3.8) 



The purpose of the transformation (3.5) is to ensure that the coefficients 6* of the partial deriva- 
tives with respect to yi in the definition (3.7) of the operator Lq are zero when i ^ d. The matrix 
d is symmetric and positive definite, but now the constant of strict ellipticity depends on b'^/b'^, 
that is, on b'^ and A, and on the constant of strict ellipticity, Aq, of the matrix a. 
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Lemma 3.3 (Local estimates for first-order derivatives of v parallel to SM). Let < r < tq, and 

let V G C2(5+) n obey 

Lqv = on , 

and assume that v satisfies 

Dv, yaD'^v e C{B:I:J and y^D'^v = on9oS+. (3.9) 
Then there is a positive constant, C = C{r,rQ,d, XQ,bo,A), such that 

hyJciB+) ^ C'll^llc(B+)' ^ d. 

Proof. We adapt the finite-difference argument employed by Brandt (1969) in [4| to prove the local 
estimates for derivatives, iJy^, when k ^ d. We let r2 := (r + ro)/2 and ra := min{(ro — r)/2, 1/2}, 
and consider the {d + l)-dimensional cylinder, 

C := {(y, yd+i) G M x M+ : y G 5+, < y^+i < rs} . 

We consider the auxiliary function, 

4>{yi Vd+i) ■■= ^ {v{y + yd+iek) - v{y - yd+iek)) , V(y, yd+i) G C, 

where C is defined above, and Ck G M'^ is the vector whose coordinates are all zero except for the 
fc-th coordinate, which is 1. We choose a constant cq > small enough, say cq = Ao/2, such that 
the differential operator, 

Lq := Lo - coyddy^y^ + c^yddy^^^y^^^, 
is elliptic on H x ]R-|_. By the definition of the function (j), we notice that 

Lg0 = on C, 

because Lqv = on . For y^ G , we consider the auxiliary function defined on C, 



V':=Ci||t;||c(^+) 



ro- 

/d-1 ^ 

yd+i{l - yd+i) + C2 Y.^y; - y^,f + yl{yd - y°) + vl+i 

\i=l / 



where the positive constants Ci, C2 will be suitably chosen below. We want to choose C2 suffi- 
ciently small that 

L^V > on C. 

By direct calculation, we obtain 

,d-l, 





= 2CiC2\\v 






= 2CiC2\\v 


)' « = 1, 2, . . . , d - 1, 


Vd+iVd+i 


= 2Ci (C2 - 


-l)ll^llc(5+)' 




= 2xdCiC2 


Mc(Bt,){^2yd-yi), 




= 2CiC2\\v 
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and so, 



= -2yrfCi||t;||^(5+) 
> -2ydCi\\v\\^^^+^ 



Co 



Co 



on C, 



using the facts that the a", for i = 1, . . . , d, and are positive constants, while 6* = 0, i 7^ d, by 
the transformation (3.5), and yd < r. We choose the constant C2 such that 

C2 <co f^a^^ + 3ra'^'^ + 2r6'^^ , 



,1=1 



so that we have 



> on C. 



Because Lqc/) = on C, the preceding inequahty yields 

Ll (±(/) - V') < on C. 

By the definition of the auxiliary function, ip, and using the fact that G and < yd+i < 1/2, 
we may choose a positive constant, Ci = Ci(r, tq, C2), large enough that 

±(/>-'0<O on^iC. (3.10) 

The portion diC of the boundary of C consists of the sets 

{Vd+i = 0, y G B+}, {yd+i = ye 5+}, and {yd+i G (0,r3), y G 

To establish inequality (3.10) along the portion {yd+i = 0} of the boundary, diC, note that (j) = 0, 
and so (3.10) holds on this portion of the boundary since "0^0. For the second portion of the 
boundary, diC, using the fact that rs < 1/2, we notice that 

yd+i(l - Vd+i) + C2 {Yl(yi - yi)^ + ydiVd - Vd) + yl+?j > rzll on {y^+i = rg, y G 5+}. 

For the third portion of the boundary, using the fact that y^ G and y G B^^ and r < r2, we 
see that on {y^+i G (0,r3), y G we have 

y,+i(l - y^+i) + C2 {^^y^ - y\? + yl(3)d - vl) + y^i j > C^^d - l)(r2 - t)\ 

Therefore, we can find a constant C3 = C3(r, rg, d) such that 

V' > C'iC3||i;||^(^+^ on {y^+i = r3, y G S+} U {y^+i G (0, r3), y G ^ifi+j 

We may choose the constant C\ = Ci{r,ro,d) large enough so that C1C3 > 1, and using the 
definition of (p, we have 

> |(/.| on {yd+i = r3, ye B+} U {yd+i G (0,r3), y G 91^+}. 

Now, inequality ( fslol ) follows. By ([3l9]) we have G C(C), and Dcj), ydD'^cf) G C(C U BqC), and 
ydD'^cj) = on SqC, where 9oC is the interior of {yd = 0} flC. Since -0 G C°°(C), we may apply the 
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comparison principle [16' Theorem 5.1] to 4> and tp on the domain C. We find that ^cj) — '0^0 
on C, and so by the definition of the function (j), we have, for ah G and yd+i £ (0, rs), 

— *— + Vd+iek) - v{yo - yd+iek)\ < Ci\\v\\^fQ+ . (1 - yd+i + C2yd+i) ■ 
^yd+i ^ '■0^ 

The preceding inequahty yields 

\vy,iy'')\<Ci\\v\\c^s+y yy'^Bt, 

for a constant Ci = Ci(r, tq, d, Aq, 6oi ^)) and this concludes the proof. □ 

Lemma 3.4 (Local estimates for higher-order derivatives of v parallel to OH). Let /c G N and 
< r < tq. T/ien t/iere is a constant C = C{k, r, ro, d, Aq, &o> such that for any v G C°°(5+) 
obeying 

Lqv = on S^, 

we /laue 

^^^1^11^(5+) < '^ll^llc(B+)' (3-11) 
for all multi-indices (3 = . . . ,/3rf_i,0) G N'^ such that |/3| < A:. 



Proof. Lemma 3.3 establishes the result when = 1. We prove the higher-order derivative 
estimates parallel to dM by induction. We assume the induction hypothesis: For any < r < ro, 
there is a constant Ci = Ci(/c — 1, r, ro, d, Aq, 6o; A), such that 

\\D^'v\\ciB+) <Cilbllc{B+)' 

for all multi-indices l3' = {(3[, . . . , /3^_i, 0) G N'^ such that \I3'\ < k - 1. Since Lqv = on 5+ , we 
also have that LqD^v = on for all multi-indices /3 with Pd = 0. We fix such a multi-index 
/3. L et A; G N be such that (3k ^ 0, and set /3' := (3 — ek- We set r2 := (r + ro)/2 and apply Lemma 
3.3 to V with < r < r2 to obtain 



for some positive constant C2 = C2{r,r2,d, Xo,bo,A). The conclusion now follows from the 
preceding estimate and the induction hypothesis applied to D^'v with < r2 < ro, since |/3'| < 
k-1. □ 



From (3.5), we have 

D^u{x) = D^v{y), y = X + ^Xd, x£ 



for all /3 G N such that 13d = 0. Therefore, Lemmas 3.3 and 3.4 give us the following estimates 
for D^u. 

Lemma 3.5 (Local estimates of higher-order derivatives of u parallel to dM). Let k £ N and 

ro > 0. Then there are positive constants, ri = ri(ro,6o,A) < ro and C = C{ro,d,k, XQ,bo,A), 
such that for any function u G C°°(-B^) solving 

Lou = onB+, (3.12) 

we have, for all /? G N'^ with (3d = and < k, 
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Proof. Let : IHI — )• EI be the affine transformation y = 4>{x) = x + £,Xd, x £ M, where ^ G M'^ 
is defined by (3.6). Let sq = so{rQ,bo, A) > be small enough such that Bf^^ C (j){B^^). Then, 
V e C°°(5+) and Lqv = on since u G C°°(S+), and Lqu = on Let si = So/2 and 
apply Lemma |3.4| to v with r replaced by si and ro replaced by sq- For any k £ N, there is a 
positive constant C = C{k, vq, d, Aq, 6o) such that for all /3 e N'' with /3rf = 0, we have 



\\dP 



(3.13) 



We now choose ri = ri(si,5o,A) small enough such that C -B^. Using the fact that 

D^u{x) = D^v{(/){x)), we obtain 

ll-^^^llc(B+ ) — ll-^'^^llc(B+ ) facts that 0(5.^) C i?J and u{x) = v{(p{x))) 

^C!MciB+,) (by ([3l3l)) 

< C||n||^^^+j (by the facts that i?^ C 0(B+) and it(x) = i;((/)(x))). 
This concludes the proof. □ 

3.3. Interior local estimates for derivatives normal to the degenerate boundary. We 

again shall use the affine transformation (3.5) of coordinates, but now with a different choice of 
the vector ^, that is 

f := -a^'^/a'^'^, Vi / d, = 0, (3.14) 
and, given a function u on H, we define the function w by 

u{x)=:w{y), y = X + £,Xd, xGM. (3.15) 



Then, by analogy with (3.7), we obtain 

Lqw := yaa'^Wy-y^ + b'wy- on 



where we notice that a*'^ = by (3.8) and the choice of t he v ector Also, we have that 
D^u{x) = D^w{y), for all /3 G N'^ with (3^ = 0. Thus, Lemma 3.5 applies to w, and we obtain a 
priori local estimates for all derivatives of w parallel to dM. 
Next, we derive an a priori local estimate for Wy^. 

Lemma 3.6 (Local estimate for Wy^). Let < r < tq. Then there is a positive constant, 
C = C(r, ro, d, Aq, &0) A), such that for any function w G C°°(i?^J obeying 



we have 



Lqw = on B^^, 



(3.16) 



Proof. Because a**^ = 0, for all i ^ d, we can rewrite the equation Lqw = on 5+ as 



VdWy^Vd + ^Wy, = f on a 



ro' 



where, for simplicity, we denote 9 := b'^ /a*^ > 0, and define / by 



d-l 



-^Wy^ onB+. 



i=l 
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We can estimate ™ terms of ) by applying Lemma 3.5 to control the supremum 

norms of Wy. and Wy.y. on B^, for all i,j 7^ d. The preceding ordinary differential equation can 
be rewritten as 

(y>y,) =yyf onB+ 
and, integrating with respect to y^, we obtain 



rvd 

ydWya{y)=l f{y\s)s^'^ ds, yeB+, 

Jo 



where denote y = (y', yd), and use the fact that 6 > 0, and Wy^ £ C{B^^^). Thus, we have 

ri/d I 

\yd'Wyaiy)\ < Wfiy', ■)\\c{[o,ya]) ^^^^ ds = -^yaWfiy', ■)\\ciio,ya])^ y^B+, 

from where it follows, by the definition of /, that 



\wy,{y)\<C \\D^^{y',-)\\cilo,yd])^ ^ ^ 



/9d=0; m<2 



for some constant C = C{XQ,bQ, A). Now applying Lemma 3.5 to estimate D^w on B^ , for all 
< r < ro, and for all /3 E N'^ with = and < 2, we obtain the supremum estimate for 
Wy^ on B^ in terms of the supremum estimate of w on B^^ . □ 

Lemma 3.7 (Local estimates for D^Dy^w with = 0). Let A; S N, and let < r < vq. Then 
there is a constant, C = C(/c, r, ro, d, Aq, 60, A), such that for any function w G C°°(5+) obeying 
(3.16) we have 

for all /? E N'^ with = and \(3\ < k. 

Proof. Since Lqw = on we also have LqD^w = on for all /3 E N*^ with jid = 0. 
Lemma 3.6 then applies with r replaced by r2 = (r + ro)/2, and gives us 



\\D^DyM\ciBt)^C^\\^^^\\c(Bt,V 



to estimate 



where Co = Co(r, ro, d, Aq, 60, A) is a positive constant. Next, we apply Lemma 3.5 
D^w and give a constant Ci = Ci(/c, r2, ro, d, Aq, &0; A) such that 

Now combining the preceding two inequalities, we obtain the a priori local estimate for D^Dy^w. 

□ 

Lemma 3.8 (Local estimate for Wy^y^)- Let A; E N and < r < ro. Then there is a positive 
constant, C = C(r, ro, d, Ao, 60, A), such that for any function w E C°°(5+) obeying (3.16) we 
have 

\\'^ydyd\\c(B+) - ^ll'"^llc{B+ )■ 

Proof. By taking another derivative with respect to yd in the equation Lqw = on we see 
that Wy^ is a solution to 



yd E d'^(^y.)y.y. + E + i^y^)y^ + + (^^Jj/. = " E 

j,j=l i=l i,j=l 



a'^Wy^yy 
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Applying the method of proof of Lemma 

d-l 



3.6 



with 6:= {b'^ + a'^'^) /a^^ and 

d-l 



Vd 



E 



a' w. 



VdViVj 



we obtain 



ydyd\\c{B+) 



< C 



E 



i=l 



\D^w\ 



c{B+y 



/3d=0,l; |/3|<3 



where C = C{d, Aq, 6o, A) is a positive constant. We can estimate the supremum norms of D^^w 
on , for ah /3 G N"^ with /3rf = 0, 1, in terms of the supremmn norm of w on B^^ with the aid 
of Lemmas 3.5 and 3.7 Now, the supremum estimate for Wy^y^ on i?^ foUows immediately. □ 



From the definition (3.15) of w, using the fact that = 0, we have 

d-l 

UxAx) = ^^kWy^{y) + Wy^{y), 

k=l 

Ux,{x) = Uy^{y), yij^d, 

d-l d-l 



[X) = ikilWy^yAv) +'^YikWy^y^{y) + 



VdVd 



(3.17) 



k,l=l 
d-l 



k=l 



U 



(x) = Y^kWy,y^^{y) + Wy^y^ (y) , \/i^d, 



k=l 



,ix) 



■w. 



for X € H. Using the preceding identities together with the estimates of Lemmas 3.6, 3.7 and 3.8 
we obtain 

Lemma 3.9 (Local estimates for second-order derivatives of u). Let vq > 0. Then there are 
positive constants, ri = ri(ro, Aq, ^O) A) < ro and C = C{ro,d,Xo,bo,A), such that for all u £ 
C°°{B+) obeying ( [3l2| ), we have 

WD^uW 



IC(S+) ^ C'll"llc{B+)' 



for all /? E N'^ with 



< 2. 



Proof. Let : IHI — t- H be the affine transformation y = 4'{x) = x + ^Xd, x G H, where ^ € M*^ 
is defined by ( 3.14[ ). Let so = ■5o('^o, Aq, A) > be small enough such that i?+ C <^(-B^). Let 
si = si{so,bo, A) < So denote the constant ri given by Lemma 3.5 applied with ro replaced by 
So- Then, the function vu defined by (3.15) has the property that w G C°°(^+) and Lqw = on 
since u E C°°(i?^) and Low = on i?^^. We apply Lemma 



3.6 



if /3 = Crf, Lemma 



3.7 



if 



f3 = Ci + ed and i ^ d, and Lemma 3^, if /3 = 2ed, to the function vu with r replaced by si and ro 
replaced by so- We apply Lemma 3.5, if /3 = or /3 = + ej, for all i,j ^ d, to the function vu 
with ri replaced by si and ro replaced by sq. Then, for any k G N, there is a positive constant, 
C = C{k,ro,d,Xo,bo,A), such that for all /3 G N"^ with < 2, we have 



(3.18) 
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We now choose ri = ri(ro, Aq, 60, A) small enough such that (j){B^_^) C Bf_^. Using (3.17), we 
obtain 

ll-^^^llc(B+ ) — ll-^^^llc(B+ ) facts that (p{B^^) C Bf_^ and u{x) = w{(p{x))) 



<C^lkllc(B+) (by (|3l8D) 

< C||n||^^^+) (by the facts that C and u{x) = w{(l){x))). 

This concludes the proof. □ 



4. Polynomial approximation and Taylor remainder estimates 

We adapt and slightly streamline the arguments of Daskalopoulos and Hamilton in [6, §1.6 & 
1.7] for their model degenerate-parabolic operators acting on functions u{t,x), {t,x) G M-(- x M?, 
to the case of our degenerate-elliptic operators acting on functions u{x), x G W^. The goal of this 
section is to derive an estimate of the remainder of the first-order Taylor polynomial of a function 
u on half-balls centered at points in dM (Corollary |4.7[ ). This result, when combined with the 
interior Schauder estimates of section ^ will lead to the full Schauder estimate for a solution on a 
half-ball centered at point in dM (Theorem |3.2[). Throughout this section, we continue to assume 



Hypothesis 3.1 and so the coefficients, a,b,c, of the operator L in (1.3) and the coefficients, a,b 



of the operator Lq in (1.4) are constant 



We let T^v denote the Taylor polynomial of degree k of a smooth function v, centered at a 
point P E M*^, and let := v — T^ denote the remainder. We then have the following analogue 
of [H Theorem 1.6.1]. 

Proposition 4.1 (Polynomial approximation). There is a positive constant C = C((i, Aq, &0) A), 
such that for any tq > 0, and any function u G C°^{B^^), there is a polynomial p of degree 1, 
such that for any r G (0, tq) we have 

„2 



\'^-P\\c{B+) 



< c 



MciB+j+^o\\Lou\\c^S+) 



(4.1) 



Proof. We first consider the case when tq = 1 and then when rg > is arbitrary. 



Step 1 (ro = 1). We let / := Lqu and we choose a smooth, non-negative, cutoff function, ^p, 
such that 



B 



1/2 



1 and \- 



m\B+ 



0. 



We fix some v > I, and let S 
solution, ui G C°°{S), to 



^ X (0,1^) as in (1.10). By Theorem B.3 there is a unique 



[LqUi 

\ui{-,i 



= 



on S, 



on 



nd-l 



Then, by setting U2 := u — ui, we see that U2 G C°°{B^^) and satisfies LqU2 = {I — tp)f on 
Notice that the definition of the functions ui and U2 differs from that of their analogues, h and 
f — h, in the proof of [6l Theorem 1.6.1]. The reason for this change is that the zeroth-order 
coefficient in the definition of Lq is zero, and so uniqueness of C°^{M) solutions to the equation 
Lqu = / on hi does not hold since we may add any constant to a solution, u. Since u = ui + U2, 
we have 



T?U2\ 



C{B+) 



< \\u2-T^U2\\c(^+^ + \\ui 



\C{B+)- 



(4.2) 
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By the Mean Value Theorem, we know that 

\\U2 - T^U2\\ci^B+) ^ Cr^\\D^U2\\c(B+y 
where C = C{d). Because LqU2 = on -6^2' apply Lemma 3.9 to U2 with r = 1/2. Then 

there are constants, ri = ri{d) and C = C((i, Aq, 6oi ^)) such that for any r G (0,ri) we have 

ll^^^2|lc-(S+) < C\\u2\\c(^B+^y 

from where it follows that 

\\U2 - Tfu2\\ci^B+) ^ Cr'^\\Mc{B+y 
Corollary A. 2 gives the estimate 

ll^l|lc(IR''-ix{0,i')) < C'||'0/||c(Kd-ix(O,i.)) < ll/llc(B+)' (4-3) 

where the second inequality follows because the support of ijj is contained in . Since u = iii+U2, 
we have 

ll'"'2|lc(B+) - ll^llc{B^) + ll^lllc(B+)' 
and so, combining the preceding two inequalities. 



Thus, we have proved that 



l"2|lc(B+) < ll'"llc(B+) + WJ IIC{S+)- 



\\u2-T^U2\\c(^B+) ^ +<^ll/llc(B+)' i^^n)- 

When r G [t'i, 1), we have, for all x G 5+, 

|u2(x) - Ti%(x)| < dll^MaWlr + \u2{x)\ + |u2(0)| 

^llc(B+) Lemma [3^ and the fact that ri < r), 

where C = C{d) \s & positive constant. Combining the cases < r < ri and ri < r < 1, we 
obtain 

||n2-r^2|lc(B+) <C^'''ll^llc(B+)+^ll/llc(5+)' Vr G (0,1), 
for a constant C = C{d, Xo,bo, A). The preceding estimate together with the identity u = ui + U2 
and (4.3) show that 

\\U - T^U2\\c^B+) ^ C (r^\\u\\c(^B+) + 11^0-^11^(5+) 

and so, the conclusion ( |4.1[ ) follows with p = T^U2, in the special case when rg = 1. 

Step 2 (Arbitrary ro > 0). When ro > is arbitrary, we use rescaling. We let u{x) := u{rox), 
for all X G Bf, and we see that {Lqu){x) = ro(Lo'u)(rox). Notice that the rescaling property 
{Lqu){x) = rQ{LQu){rQx) does not hold in this form if the zeroth-order coefficient of Lq is non-zero. 

We apply the preceding step to u with r replaced by r/rg. Then, there is a polynomial p such 
that 



which is equivalent to 

ll^-^'llc(B+) ^ (^lkllc(B+) +ro||LoM||^(^+)^ , 
where we set p{x) := p{x/ro). We notice that the polynomial p depends on ro, but not on r. 
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The proof of Proposition 4.1 is now complete. 



□ 



Proposition 4.1 is used to obtain the following analogue of fGl, Theorem 1.7.1]. 

Proposition 4.2. For any a G (0, 1), there is a positive constant S = S{d,XQ,bQ,A), such that 
for any u G C°°{Bf) with T^u = 0, we have 



sup 

o<r<i r 



C(B+) ^ ^ 



l+a 



\C(B+) 



+ ^ + sup 



\\Lou\ 



C{B+) 



0<r<l 



(4.4) 



Proof. Because T^u = and u G C°°{B^), it follows that the quantity on the left-hand side of 
the inequality (4.4) is finite. In addition, T^u = implies that Lqu{G) = 0, and so we also have 



ll^oull 



Let r=i, G (0, 1] be such that 



and we define for simplicity 



sup 

0<r<l 



\U\\ 



\C{B+) 
sup ^-p — - 

0<r<l ' 



< oo. 



Q ■= \W\ 



C{B+) 



+ sup 

0<r<l 



ll^onll 



C{B+) 



(4.5) 



We let S (depending on u) be such that 



SQ. 



(4.6) 



It is sufficient to find an upper bound on S, independent of u, to give the conclusion (4.4). 

Let q and s be positive constants such that 0<g<r*<s<l. We apply Proposition |4. 1| to 
u with r replaced by q and r* and rg replaced by s. Then, we can find a degree-one polynomial, 
p, such that 

,2 



\U — p\ 



\C{B+) 
C{B+) 



< C 



Mc{B+)+4Lou\\c^g, 



\u — p\\^,D+ ^ < C 



U 



\c{B+) +4Lou\\c{B^ 



But 



\\P\ 



CiB+) 



<C-\\p\ 



c{B+y 



for some positive constant C, depending only on d. We can then estimate 

,, ,, 

< C- 



^C{B+) ^ ^ „ V''" Ph{B+) + \MciB+) 



and using (4.7) and the fact that q < r^:, we obtain 

^2 



\\P\ 



C{B+) 



\c{B+) + 4^0'^\\c{B+) + ll^llc(B+) 



From 



IC(B+)' 



(4.7) 
(4.8) 



(4.9) 
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and (4.8) and (4.9), we see that 

^2 



\u\ 



CiB+) 



< c 



< C 



2 ll'"llc(B 



+ ) + 5|l^on|lc(i?+) + -^MciBt) + — ll^0^^llc(B+) + -\Mc(B+) 



Tl II II r^s 



+ -^||Lom||^(^+) 



where we have used the fact that q < < s to obtain the last inequahty. We divide by r;[+° and 
find that 



11^(5+ ) 

' * 



< c 



C(Bt) 



ll"IIC{B+) S ( S ^ 



+ 



\\Lou\ 



C(B+) 



From the preceding inequality and definitions (4.5) of Q and (4.6) of S, we deduce that 

' r \ 1— Q 



SQ<C 



+ 



SQ + C 



s s 



q \r 



Q, 



By choosing r^:/s and q/r^: small enough, we obtain a bound on 5 depending only on C 
C{d, XQ,bo, A). Hence, the estimate (4.4) now follows. 



□ 



We apply Proposition |4.2| to RiU := u — T^u. Note that LqT^u = {Lqu){0) and so 

Lo {u - T^u) = Lou - (Lou)(0) = ii[|Lou, 



because XdD^u = on dM and the zeroth-order coefficient of Lq is zero. Thus, Proposition 4.2 
yields the following analogue of [6, Corollary 1.7.2]. 

Corollary 4.3. For any a G (0, 1), there is a positive constant S = S{d,XQ,bQ,A), such that for 
any u € C°°{B^) we have 



ll^l«llc(B+ 



' ^ ^ P?«llc7(B+) + sup 

0<r<l V 0<r<l 



||i?[]Lon|| 



C{B+) 



Using the inequality (2.3), 

< 2s(x,0)^ VxGM, 

where we recall that the cycloidal distance function, for all G H, is given by (2.1), 

we see that there is a positive constant, C = C{a,d), such that 



II^O^O«|lc(S+) ^ ^ 

sup < C sup 

0<r<l 0<r<l 



Lqu{x) — Lqu{0) 



s2"(x,0) 

< C [Lou](j2c,(^§+) ■ 



C(B+) 



Therefore, Corollary 4.3 gives us the following partial analogue of [6l Corollary 1.7.5]. 

Corollary 4.4. For any a G (0,1), there is a positive constant C = C{a,d, Xo,bo,A), such that 
for any u G C°°(i?^) and < r < 1, we have 



Next, we improve the estimate in Corollary 4.4 with the following analogue of [HI Theorem 
1.7.3 & 1.7.6]. 
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Proposition 4.5. For any a E (0,1), there is a positive constant C = C{a,d, XQ,bo,A), such 
that for any u G C°°{B^) we have 



\T?u\ 



C{B^) 



+ ^ < C 



\c{B+) + [^0?^]c-{B+ 



Proof. Because is a degree-one polynomial, there is a positive constant C = C{d) such that 



1 "IIC(B+) 



< 



1 "IIC(B+)' 



Vr G (0,1]. 



By Corollary 4.4, we have for all r G (0, 1], 



\R>\\ciBt) < Cr'^""^' (hi 



By combining the preceding two inequalities, we find that 



C{B+) 







< 


-( 




r V 






< 


(- 




V r 



\c{B+) + ll-^i^llc(B+; 



+ 



C(B+) 



+ C7r°/2 II Tf- 



lc;(B+) 



By choosing r small enoug h so that Cr"/2 < 1/2, 

we obtain the conclusion. 



□ 



Proposition 4.5 implies the following special case (r = 1) of [6, Corollary 1.7.8]. 



Corollary 4.6. For any a £ (0,1), there is a positive constant C = C(a, d, Aq, &0) ^)! such that 
for any u G C°°{Bf) we have 



\RU\ciBt) < C 



u 



Corollaries 4.4 and 4.6 yield the following analogue of [6, Corollary 1.7.8]. 

Corollary 4.7. For any a G (0,1), there is a positive constant C = C{a, d, Xo,bo, A), such that 
for any u G C°°{B^) and < r < 1, we have 



C(B+) 



< Cr^+^l^ 



\C{B+) 



+ [Lqu] 



Cf{B+) 



5. SCHAUDER ESTIMATES AWAY FROM THE DEGENERATE BOUNDARY 

In this section, we use a scaling argument to obtain elliptic Schauder estimates away from 
the degenerate boundary analogous to the parabolic versions of those estimates in [Gj §1.8]. Our 
argument is shorter because we only aim to obtain the estimates in Lemma 5.1 and Corollary 



5.3 Even though these estimates are weaker than their analogues |6], Corollary 1.8.7] and [U 



Corollary 1.8.8], respectively, they are sufficient to obtain the full Schauder estimate (3.4) in 
Theorem |3.2[ The estimate (3.4) is proved using a combination of the Schauder estimate on balls 
Br{xo) (s in which we prove in this section, and the results of ^ The proof of Proposition 6.1 



uses Corollary 5.3 which is derived from Lemma 5.1 We have encountered a similar situation 
in the proof of Holder continuity along dM. of a weak solution to the Heston elliptic equation in 
|17[ Theorem 1.11]. Throughout this section, we continue to assume Hypothesis 3.1 and so the 
coefficients, a,b,c, of the operator L in (1.3) and the coefficients, a,b, of the operator Lq in ( |1.4[ ) 
are constant. 

For any r > 0, we let Qr denote the point re^ G H. We have the following analogue of [HI 
Corollary 1.8.7]. 
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Lemma 5.1. For any a G (0, 1) and positive constants /i and A such that < /x < A < 1, there 
is a positive constant C = C{a, fi, X,d, Xq, A), such that the following holds. For any function 
u G C°°{Bxr{Qr)), we have 



CS'iB^riQr)) 



+ - 



1 



.a/2 



II^O^II 



C(BAr(Qr)) 



(BxriQr)) 



(5.1) 



Remark 5.2. The estimate in [H Corollary 1.8.7] does not contain the term ||Lo?/||c(^^^(q^)) 
appearing on the right-hand side of our interior estimate (5.1 ). However, our estimate is sufficient 
to give the Schauder estimate (3.4) in our Theorem |3.2[ 



Proof The result follows by rescaling. We denote x = {ry' ,r + ry^) G H, where we recall that we 
denote y = (y', y^) G H = x M+, and define 

v{y) = u{x), \/y e Bx- 

By the hypothesis u G C°^{Bxr{Qr)), it follows that v G C°°{Bx) and v is & solution to the 
strictly elliptic equation 

+ VvyM = rf{y), yy = {y',yd) g Bx, 

where f{y) := f{ry',r + ry^), for all y G Bx, and / := Lqu. By the interior Schauder estimates 
[301 Theorem 7.1.1], there is a constant C = C{a, fi, A, d, Aq, A), such that 

(5.2) 



(5.3) 



D vWc^^B^^ < C (\\v\\c(B^)+r\\f\\c<.^B^) 



By direct calculation, we obtain 



\M\C{B^) - iMciBUQr))' 

\\f\\c{B^) = ll/llc(BAr(Qr))' 

[f]c-(B.) < Cr^'^ [/] 



Cf (BAr(Qr)) ' 



where C = C{a). To see the last inequality, recall that x = {ry',r + ry^), for all {y',yd) G Bx- 
For any G Bx, for i = 1,2, we have 



\f{y')-f{y')\ \f{x')-f{x')\six\x^r 



\yl _ y2|a 



s{x^,x'^)°' \y^ — 2/^1° 



By (2.1), we notice that 

s{x^,x'^) 



I 1 2 1 

r\y - y \ 



\y'-y'\ ./^(2 + yi + y2 + |yi_y2nir-r 



r 
2' 



and so, by letting C = 2 we obtain 
We also have 



^Cf{B^.r{Qr)) 



<^^-(l+a/2)||^2^ 



(5.4) 
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for a constant C = C{a,d). To establish ( |5.4[ ), we only need to consider quotients of the form 

where x^, G B^r, and all their coordinates coincide, except for the i-th one, where i = 1, . . . ,d. 
We only consider the case when i = d, as all the other cases, i = 1, . . . , d — 1, follow in the same 
way. Recall that we denote x = {ry' , r + ryd), for all y £ B^. We obtain 



x\D^u{x^) — x^D^u{x 
s{x^,x'^)'^ 



< 



\D\(x^)\ + x'^. 



\D^uix^) - D'^u{x^)\ 



Using the definition of the cycloidal metric (2.1), and the fact that D'^u{x) = l/r'^D'^v{y), for all 
X £ Bfj_r, we see that 



x^L>^m(x^) - x^L>^u(x^)| 
s(xi,x2)° 



< (24 + 2x^) 



2^°/2 1^1 _ ^2|1-qJ_ 



+ X 



,1 \DMy')-DMy^)\\y'-y 



2la 



yl _ y2\a s(x^, X^)' 

2„ 



|yl_y2|a 



where we used the fact that x^ < r, for all X"'^,x^ G B^,-. We also have by (2.1) 



which implies that 



\y'-y'\ 
s(x^, x^) \x\ — x^l 



-l ZjMxl)-xlZ)Mx^)| ^ ^^-(l+./2)||^2 



X 



\^xlMx\-xl\<Cr-^l\ 



S(xl,x2)" 



IC?(Bm)' 



for a constant C = C(a). Now, the inequality (5.4) follows immediately. 
Using the preceding estimates (5.3) and (5.4), it follows by (5.2) that 



where we substituted Lqu for /. 

To obtain the estimate for the Holder seminorm of Du, we proceed by analogy with the 
argument for x^D'^u. □ 

We have the following analogue of P, Corollary 1.8.8]. 

Corollary 5.3. For any a £ (0, 1) and positive constants n and X such that < /i < A < 1, there 
is a constant C = C{a, fi, X,d, Xq, A), such that for any function u G C°^{Br{Qr)) we have 



|i?rn| 



Proof. As in case of the inequality preceding |6l Corollary 1.8.8], we have 



Cf{B^r(Qr)) 



for a constant C = C{d). Thus, the conclusion follows from Lemma 5.1 and the preceding 
inequality. □ 
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6. SCHAUDER ESTIMATES NEAR THE DEGENERATE BOUNDARY 

In this section, we use the results of the previous sections to prove our main a priori interior 
local Schauder estimate (Theorem 3.2) for the operator L on half-balls centered at points in the 
"degenerate boundary", dM.. Throughout this section, we continue to assume Hypothesis 3.1, and 
so the coefficients, a,b,c, of the operator L in (1.3) and the coefficients, a,b, of the operator Lo 
in (1.4) are constant. 

We begin with an analogue of [6^, Theorem 1.9.1]. 

Proposition 6.1. For any a G (0,1), there is a constant C = C(a, d, Aq, 6o, A), such that the 
following holds. For any function u G C°°{B^) and any < r < 1/2, we have 



Proof. We choose = 1/4 and A = 1/2 in Corollary 5.3, We consider the points Qr '■= re^ and 
P := O £ W^. Let p := T^''^ — u, where we recall that T^' ti is the second-degree Taylor 
polynomial of u at Qr, and Tf u is the first-degree Taylor polynomial of u at P. Then, we also 
have that p := Rfu — F^'^u, where we recall that E^^'u is the remainder of the second-degree 
Taylor polynomial of u at Qr, and R^u is the remainder of the first-degree Taylor polynomial of 
u at P. There is a positive constant, C = C{d, fi), such that 

C 

\D^p\ < ^\\p\\c{B,AQr))^ 
which implies, from the definition of p, that 

C 



(Qr))- 



(6.1) 
(6.2) 



Corollary 4.7 applied to Rfu gives 

and the interior Schauder estimate in Corollary 5.3 applied to Rfu yields 

+ r\\LoR^u\\c^B,AQr))+'^'^''^^ t^o^^""] C?(B..(Q.)) 

for a constant C = C{a, /x. A, d, Aq, A). We notice that LoRfu = Lqu — {Lou){P), from where it 
follows that 



using the fact (2.2) that s{x^,x'^) < Ix"*^ — for all G H, and also that 



[LoR( 



u 



[Lqu] 



(6.3) 
(6.4) 



The preceding three inequalities, together with (6.2), give us 



\RfRlU\\c{B,,{Q,.)) < Cr^^''/^ {\\u\\ciB+) + [Lou]c?{Bt 



where we used the fact that B^r{Qr) C when < r < 1, for all < // < 1. Notice that 
R2'^u = i?^''i?f n, and so the preceding estimate becomes 



\R'. 



■-7,11 - ^ r<r'^+a/2 I I 



C{B+) 



[Lqu] 



C^{B+) 
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The conclusion now follows from the preceding estimate, and inequalities (6.2) and (6.1). □ 

Note that the definition ( |2.1[ ) of the cycloidal metric gives 

s{{x,Xd),ix',0)) = y/xd/2, y{x',Xd)£M, 

and hence, via Proposition 6.1, we obtain the following analogues of Theorems 1.9.3 & 1.9.4]. 

Corollary 6.2. For any a £ (0,1), there is a constant C = C{a,d,Xo,bo,A), such that for all 
<Xd< 1/2 and x' G R'^~'^ , and any function u G C°° {Bf {x' , 0)), we have 



\xdD'^u{x',Xd)\ < Cs {{x',Xd), {x',0))°' (\\u 



\C{B+{x',0)) 



Corollary 6.3. For any a G (0,1), there is a constant C = C{a,d,Xo,bo,A), such that for all 
<Xd< 1/2 and x' £ R'^^^ , and any function u G C°° {Bf {x' , 0)), we have 



\Du{x',Xd) - Du{x',0)\ < Cs {{x',Xd), (x',0))" (\\u\\c{B+{x',o)) + i^o'^]ci{B+{x' 



0)) 



Proof. Following the proof of [HI Theorem 1.9.4], using Proposition 6.1 and translation-invariance 
with respect to x' G R'^^'^ to obtain the second inequality, we have 



\Du{x' ,Xd) — Du{x' , 



< 



Du^^ix',t)\dt 



<C \\u 



\ f^"^ /2 1 
lc(B+(x',0)) + [-^O^lcf (B+{x',0))j 



C \\u\\ 



IC(B+(x',0)) + [-^O^^Cf (B+(x',0)) J -^d 

Using the fact that s{{x',Xd), {x',0)) = y^Xrf/2, we obtain the conclusion. 



a/2 



□ 



Next, we use Lemma 5.1 (for estimates away from dM) and the Taylor remainder estimates in 
Corollary 4.7 (for estimates near dM) to prove the following analogue of [6l Theorems 1.9.5]. 

Proposition 6.4. Let a G (0, 1), and < r < 1/4, and < < 1. Then there is a constant, 
C = C(a, ;U, d, Ao, A), such that for any function u G C°°{B^), we have 



Proof. For simplicity, we denote v := R^u = u — Tf u. We notice that 



\Du] 



C?(5^r(Qr)) 



\Dv 



C?(B^,r{Qr)) 



and [xdD'^u]c^^B,AQr)) = I^d^^^^lcf 



))' 



and hence we only need to estimate [Dv]fja/^ 



(BMr(Qr)) 



and [xdD'^v](^a(^^^^(^Q^'^y The proof is similar 
to the proof of Proposition 6.1 The interior Schauder estimates in Lemma 5.1 applied to v with 
A = (1 + ^)/2 yield 

1 „ „ 



[Dv_ 



Cf{B,,r{Qr)) 



lCf(B„.(Q.)) 



< c 



\v\ 



C(BAr(Qr)) 



+ -n^\\LM\c{B^,.{Qr)) 



+ [Lqv] 



Cf(BAr(Qr)) 



j.l+a/2 I 
1 

for some constant C = C(a, /x, d, Aq, A). The conclusion now follows from the preceding estimate 
and inequalities (6.2) applied on Bx^iQr) instead of B^ (notice that Bxr{Qr) C By^, since 
< r < 1/4), together with and ( [6^ . □ 

Next, we have the following analogue of [H Theorems 1.9.7 & 1.9.8]. 



REGULARITY OF SOLUTIONS TO DEGENERATE-ELLIPTIC EQUATIONS 



29 



Proposition 6.5. For a G (0, 1), there are constants 7 = 7((i) G (0, 1), and C = C{a, d, Aq, 601 
such that for any function u G C°°{Bf), we have 



[Du] 



Cf(B+) 



+ [xdD\]^^^^+^ < C 



\u\ 



C{B+) 



+ [Lqu] 



Cf{B+) 



(6.5) 



Proof. We combine the arguments of the proofs of [U Theorems 1.9.7 & 1.9.8]. Let G , for 
i = 1,2, where 7 wiU be fixed below. We may assume without loss of generality that x"^ > x^. 
We consider two cases. 

Case 1 (x^ and close together relative to dist(-, 9IHI)). If |xi — xi| < x^/4, then x^ G B^ii^{x'^), 
and the estimate (6.5) follows if we assume < 7 < 1/2 and apply Proposition 6.4 with /i = 1/4 
and r = x^. 

Case 2 (x^ and x^ farther apart relative to dist(-, 9IHI)). We next consider the case when 

|x^ -x^l > x^/4. (6.6) 

Writing x = (x, x^) G 



X ]R_|_, we define the points, 

x3:=(x\0) and x^:=(x^O), 



X 



5 ._ 



(x"*^, r) and x^ : = 



X , r 



'I h 



where the positive constant r will be chosen below. Notice that when (6.6) holds, we have 

s(x"'", x^) > -a/^!^, 



by the definition (|2.1|) of the cycloidal distance function. By the definition of the points x*, for 

(6.7) 



i = 3, 4, and the fact that s((x', x^), (x', 0)) = Y^Xrf/2, we see that 

six^.x^) > 8s(x\x^), 
s(x^,x^) > 8s(x^,x^) (since x^^ > x^). 
Let f denote Du or XdD'^u, and consider the difference 

u(x^) - v(x^) = (u(x^) - f (x^)) + (t;(x^) - v{x^)) + {v{x^) - v{x^)) 
+ {v{x^) - v{x^)) + {v{x*) - v{x^)) . 



(6. 



Using the distance inequalities (6.7), we find that 

< 8° 



|t'(x^) — f (x'^)| 

S(xl,x2)" 

|t>(x^) — f (x^)| 
s{x^,x'^)°' 



< 



|f (x-*-) — f (x^)| 

S(xl,x3)" 

|f (x^) — ^(x^)! 
s(x^, x^)" 



By Corollary 6.2, if v = XdD u, and Corollary 6.3, if v = Du, we obtain 
i'(x^) — ?;(x^)| |?;(x^) — f (x^) 



+ 



S(xl,x2)" 



(6.9) 



for a constant C = C(a, d, Aq, bo; A). 

We now let r := i?s^(x^,x^), where the constant B will be chosen below. Using the fact that 
s(x^, x^) = \/r/2 and definition of x*, for z = 3, 5, we obtain 



|u(x^) — f (x^)| 



s(x^, x^)" 



{B/2y 



y/2 I^(X^) -^(X^)| 



s(x^, X 



5^a 
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Because G B^, for i = 1,2, and due the inequality (2.2), we can choose the constant B :- 
1/(47) such that 

r = Bs^{x^,x'^) < B\x^ - x^l < B-f < 1/4. 
We apply Corollary 6.2, when v = x^iD'^u, and Corollary 6.3, when v = Du, to obtain 
If (x'^) — v{x^) I 



The inequality, 



<CiB/2r/'[\\u\\^,^^+^ + [Lou 



\Cf{B+) ■ 



^-2^ <C{B/2T'^ ( ll^llc(5+) + [^0^1^.(54 



s{x^,x'^) 



(6.10) 



(6.11) 



follows by the same argument used to obtain the estimate (6.10). 
Using (6.6) and the assumption x\ > x^, we see that 

3 



\x^ — x^l < Ix"'" — x^l < 



x^ — X^P 



2 x^ + x^ + — x^l 2 



s\x\x^). 



Recalling that B = 1/(47) ^iid ^ = i?s (x , x ), we have 



Ix"^ — x^l < -—Bs^ix^ , x^) < 67r. 



Next, we choose 7 = 1/24, and so 



Ix"*" — x'^l < r/4. 



for all X* = (x^, • • • , x^) G B^, for i = 1, 2. Because |x^ — x^| < r/4, we may apply Proposition 



6.4, with 11= 1/4, to obtain 

\v{x^) - v{x^)\ ( 

Again using the definition r := i?s^(x^,x^), we notice that 

S(X^X6) < "f'l < '^S(x\x2), 

V 2r 4 



and so the preceding two inequalities yield 

It'(x^) — ?;(x^)| 



S^i^X^^x"^) - ^ \^^"'^^C{B7) ' L^-u-JC'^fBTW • (6.12) 

Combining the estimates (6.9), ( |6.10 ), ( |6.11 ) and (6.12) gives us the estimate (6.5), when condition 
( [6^ holds. 

The conclusion now follows from the two cases we considered. □ 

By analogy with [6, Corollary 1.9.9], we have 

Proposition 6.6. For any a G (0,1), there are positive constants 7 = ^{d) G (0,1), and C = 
C{a,d, Xo,bo, A), such that the following holds. If u G C°°{Bf), then 



(6.13) 
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Proof. Le t 7 = 'y{d) £ (0 , 1) be as in Proposition 
Corollary 



6.2 



6.5 



The bound on XdD^u follows from 
Proposition 6.5 gives us the estimate (6.13) for the Cg{B^) Holder seminorms 
of Du and XdD'^u. We only need to show the bound on of Du, namely that there is a constant 
C = C{a, d, Ao, 60, A), such that 



C{B+) 



< C lln 



\C(B 



+) + [^ou]c.„(B+ 



We follow the argument of ||6i p. 932]. Let € be such that \Du{x^ 
by Proposition 



\Du\ 



6.5 



C{B+) 



(6.14) 
. Then 



we have, for all x £ B2. 



\Du{x) - Du{x^)\ < Co Mhllc(B+) + [Lm]cf{Bi 



for a constant Co = Co(a, d, Ao, ^o. A). 

Let > 2 be a positive integer such that 



\D^\\c(B+ 



) > NCq + [io'u]c-(B+)) ■ 



Estimate (6.14) will follow if we can find an upper bound on iV, independent of u. The preceding 
two inequalities give 

\Du{x)\ > (iV-l)Co(||n||^(5H.) 

and the Mean Value Theorem yields 



\u{x) - u{x^)\ >\x- x^\{N - l)Co ( \\u 



\C(B+) 



[Lqu] 



C?{B^) J ' ^^^^ ^ -^7 



Choosing x € B^ such that |x — x^\ > 7/2, we obtain a contradiction with (6.5 ) if is too large. 
Thus, dell) fohows. □ 



We have the following corollary of Proposition 6.6 



Corollary 6.7. For any a G (0,1), there are positive constants, 7 = "y{d) G (0,1) and C 
C(a, d, Ao, 60, A), such that for any r > the following holds. If u £ C°°{B^), then 



I^llc2+"(B+,) 



(6.15) 



Proof. Let 7 = jid) G (0, 1) be as in Proposition 6.5 We set v{x) := u{rx), for all x G B^ . The 
estimates of Proposition |6.6| applied to v give us 



where C = C(a, r, d, Ao, ^0, A). The dependency of the constant C on r follows as in the proof of 
Lemma 5.1, and so we obtain (6.15). □ 



We now generalize Corollary 6.7 to allow for any 7 G (0, 1) and make explicit the dependency 
of the constant C appearing in (6.15) on r and 7. 

Corollary 6.8. //a G (0,1), then there are positive constants, p = p{a) andC = C{a,d, Xo,bo,A), 
such that, for any r > and 7 G (0, 1), the following holds. If u £ C°°{B^), then 



< C((l - 7)^-) " {\Mc{B+) + [^0'"]c-(B+)) ■ 



(6.16) 



32 



P. FEEHAN AND C. POP 



Proof. Let £ (0, 1) be the constant appearing in the hypotheses of Corollary 6.7 We consider 
two cases, < 7 < 7^ and 7 > 7^, and clearly we only need to consider the second case. Our 
proof of (6.16) now follows standard covering argument. Let t := (1 — 7)r/2 and divide the 
half-ball, B^^., into the two regions. 



u, ■= B± n 



7r 



x(0,7dV2) and U2:=B+\Ui 



We cover Ui by a finite number of half-balls, B^^^{x'^), centered at points G doB^{x^) and we 



apply the estimate (|6.15|) to obtain 

\c, 



< (ii^i 



C{B+) 



+ [Lou] 



Cf{B+) 



<C((l-7)r)-(i+"/2) [\ 



\C{B+) 



+ [Lqu] 



C?{B+) 



(using t = {1 — 7)r/2) 



In the region U2, the operator Lq is strictly elliptic (because Xd > 7di/2 > 0), and so we may 
apply |26| Corollary 6.3]. Using the elliptic analogue of the parabolic estimate [18] Proposition 
3.13], there is a positive constant p = p{a) such that 



)<C((l-7)r)-p(||n||^(^+) + [Lon 



JCf{B+) 



□ 



Estimate (6.16) follows by combining the preceding two inequalities. 

Proof of Theor em\3.S\ We combine the localization procedure in the proof of \30\ Theorem 8.11.1] 
with Corollary 6.8 We divide the proof in two steps. Set R := (r + ro)/2. 

Step 1 (A priori estimate for u G C°^(i?,^)). Consider the sequence of radii, {rn}n>i C [r,R), 
defined by ri := r and 



r+{R 



n—l ^ 

y-- 

k=l 



Vn > 2. 



(6.17) 



Denote Bn := (x'^), for all n > 1. Let {ipn}n>i be a sequence of Cg°(]HI) cutoff functions such 
that, for all n > 1, we have 1 with = 1 on Bn and (fn = outside Bn+i- Let 

Vn > 1. (6.18) 



a. 



\Ufn\\r^+'^ 



By applying the estimate (6.16) to ucpn with r = r„+i and 7 = rn/rn+i, we obtain 
an < C(r„+i - r„)-P [\\u^Pn\\c(Br^+l) + I^o^^lcf 

< C(i? - r)-P2("-l)P (ll^V'nllcCB^+i) + ll^^^llci'CB^+i) + ll'"'/'n+l|lc|"{B„+i) } ' 



where the last inequality follows from the fact that L = Lq -|- c by (1.3) and (1.4) and employing 
(6.17). The interpolation inequalities (Lemma C.2) give, for any e > 0, 



hfn+l\\cS'{B„+i) < ^(^n+l + Ce '^\\uipn+l\\c{B,,+^) (by ( |6.18D ), 

where C = C{d,a,R) and m = m{d,a) are positive constants independent of e. Choosing 
e := 5C-i2-("-i)P(i? - r)P, we obtain, for all S > 0, 

an < San+i + c(^iR- r)-P2("-i)P + 5-™(ii - r)-p(™+i)2("-i)(™+i)pj 



\Lu\ 



Cf{B+{x^)) 



+ \\u^ri+l\ 



C?(B+(x<^))) ' 



and now the estimate (3.4) follows as in the proofs of |30| Theorem 8.11.1] or |181 Theorem 3.8]. 
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Step 2 (A priori estimate for u G C^+"(S+ Clioose a sequence {un}nm C C°°(i?^(x°)) 
such that tin —5- in 0'^'^°' {B^{x^)) as n — )■ oo. Applying the estimate (3.4) to each u„ and then 
taking the limit as n — t- oo, yields the a priori estimate (3.4) for u G C'j'^"' {B^^{x^)) . 



This concludes the proof of Theorem 3.2 



□ 



7. Higher-order a priori Schauder estimates for operators with constant 

coefficients 



In this section, we prove a higher-order version of Theorem 3.2 our basic a priori local interior 
Schauder estimate, and a global a priori global Schauder estimate on a strip (Corollary 7.2), both 



when L has constant coefficients. Throughout this section, we continue to assume Hypothesis 3.1 



and so the coefficients, a,b,c, of the operator L in (1.3) and the coefficients, a,b, of the operator 
Lq in (1.4) are constant. 

Theorem 7.1 (Higher-order a priori local interior Schauder estimate when L has constant coef- 
ficients). Assume the hypotheses of Theorem 3.2 and let A: G N. If u G Cs'''^~^'^{B^{x^)), then 



\u\\^k,2+afp,+ -^ < C ( ||Lu||^fe,c«|,p+ ^^o^^ + ll^l 



(7.1) 



where C now also depends on k. 

Proof. Choose ri := (r + ro)/2 G (r,ro). For any multi-index /3 G N'^ with |/3| := /3i + - • ■+/3^ < k, 
direct calculation yields 

where we write f3o := /3 — f^dCd and, for / G N, 



V G C~(M), 



(7.2) 



id 



Vx, + [h'' + la'^'^ ) + cv. 



Note that L(o) = L. To prove (7.1), we see by Definition 2.3 that it suffices to establish 



\D"u\ 



< C (ll^^^llcJ,"(B+{xO)) + 11^ 



(7.3) 



for any multi-index /3 G N"' with |/3| < A;, where C has the dependencies given in our hypotheses. 

Theorem 3.2 yields (7.1) when k = 0. Therefore, as an induction hypothesis for k, we assume 
that (7.1) holds with k replaced by any Z G N in the range < / < A; — 1 and we seek to prove 
(7.3) and hence (7.1) by induction on I when |/3| = k. 



e case 


Pd = 0, 


SO LDl^v 


- ^ ( 


WLD^u 


lcf(B+(xO 


- ^ ( 


WDI^Lu 


lcf(B+(xO 


- ^ ( 


\\Lu\\„k 




- ^ ( 




-(B+(xO)) 


- ^ ( 







+ \\D^u\ 
+ \\D^u\ 



C(5+{xO)) 



(by (3.4) 



+ ll^llc,*-°(B+{xO)) 



(by (jT^) 
(by Definition [2^ 



(by Definition 2.3) 
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where the final inequality follows by induction on / and the a priori Schauder estimate (7.1 ) with 



k replaced hy I = k — 1 (and r replaced by ri). Since ri < ro, we can combine terms and obtain 



(7.3) in the case /3d = 0. 

Now we consider the case < f3d < k and argue by induction on j3d- As an induction hypothesis 



for f^d, we assume that (7.3) holds when < (3d < k — 1. For (3d in the range I < (3d < k (and 



thus \(3o\ < k — 1), we have 

ll^^"llc?+"(B+) 

< C {WhMD^lcHB^.iccO)) + P^«llc(B+ (.0))) (by 0) 



»^jllCf(B+{xO)) + 11-^ ""I 



CiB+ixO)) 



2.3 



where the penultimate inequality follows from (7.2) and the final inequality by D efi nitio n 2.2 of 
our Holder norms. Because C2+°(B+(x°)) ^ aP(-B+(x°)) by Definitions [^i] and 
that 



we see 



(by induction on (3d and (7.3) since /3d — 1 < A; — 1 



< C ( (^0)) + ll«llc,*>°{B+ (xO)) + ll^llc(B+ (xO)) 



(since ri < tq) 



Cf'"(B,+ (x")) ^ ll"'llc'=-i.2+"(B+(xO)) + ll^llc(S+{xO)) 



< C ( (xO)) + ll^^llc*-''"(B+ (xO)) + ll^llc(B+ (xO))y ' 

fc-l,2+a 



where the penultimate inequality follows from the embedding Cg ~ ' "(-B^ (x^)) ^ Cg '"(-B^ (x^)) 
implied by Definitions 2.2 and 2.3 and the final inequality follows by induction on / and the a 
priori Schauder estimate (7.1) with k replaced by / = A; — 1 (and r replaced by ri). Again, since 
ri < ro, we can combine terms and obtain (7.3) in this case too. □ 



Let ly > and let S = M.'^^'^ x (0, u), as in (1.10). Theorem 7.1 together with a priori estimates 
for strictly elliptic operators in [26^ §6] now imply the following global Schauder estimate on 
strips. 

Corollary 7.2 (A priori global Schauder estimate on a strip when L has constant coefficients). 

For any a £ (0, 1), constant z/ > 0, and k £ N, there is a positive constant, C = C{k, a, v, d, Xq, bo, A), 
such that the following holds. If u £ Cs''^~^°'{S) and u = on diS, then 



and, when c > 0, 



(7.4) 
(7.5) 
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t>d-l 



(0,r/4)c [jB^^^i 



neN 



Using the a priori interior local Schauder estimate (3.4) on each half-ball B^{x"'), we obtain 

+ \\u\ 



\\U\\„k,2 + a,^ + 
s ^ r / 

By applying a standard covering argument to the strip Sq : 



C(Bt{x")), 

X (0,r), we find that 



By \26\ Lemma 6.5 &: Problem 6.2] and a similar covering argument, there is a constant S > 
such that, if Si := x (i/ — (5, z^), we have 



Setting 52 := R'^'^ x (r/4,z^ - 5/2) and now applying [26, Corollary 6.3 & Problem 6.1] and a 
covering argument, we obtain 



By combining the preceding three estimates, we obtain (7.4) and by appealing to Corollary A. 2 



we obtain (7.5) 



□ 



8. A PRIORI Schauder estimates, global existence, and regularity for operators 

WITH VARIABLE COEFFICIENTS 



In { 8.1 we relax the condition in Hypothesis 3.1 that the coefficients, a, b, c, of the operator L in 



(1.3) are constant, which we assumed in sections [3|[4|[5| and [6| to prove a generalization (Theorem 
8.1 ) of our Cg^" a priori Schauder estimate (Theorem 3.2 ) from the case of constant coefficients. 



a,b,c, to the case of variable coefficients. We then prove Theorem 8.3, extending the preceding 
C^^" a priori Schauder estimate to a Cs a priori Schauder estimate for arbitrary A; G N. This 
allows us to complete the proofs of Theorem 1 1 . 1 1 and Corollary 1.3 In § refsubsec: Regularity, we 



prove our global Cs' "(5) existence result on strips, S, and hence a Cg' "(^^(x''))-regularity 
result, Theorem 8.4l on half-balls, (x^). We conclude the section with the proofs of Theorems 



1.8 and 1.11 , and Corollary |1.13 



^.1. A priori Schauder estimates for operators with variable coefficients. We begin 



with a generalization of Theorem 3.2 to the case of variable coefficients. 



Theorem 8.1 (A priori interior local Schauder estimate when L has variable coefficients). Let 

a € (0, 1) and let tq, Aq, bo, A be positive constants. Suppose that the coefficients a^^ , ¥ , and c 
of L in (1.3) belong to C^{B^ (x^)), where x^ £ dM, and obey 



Mci{B+ixO)) + Mci(B+(xO)) + Mci{B+{xO)) < ^' 

b'^ybo ondoB+ix"^), 



(8.1) 
(8.2) 
(8.3) 
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Then, for all r G (0, tq), there is a positive constant C = C{a, r, ro, d, Xq, bo, A) such that, for any 
functioi^u G C^^"(i?^(x'')), we have 



lc?+"(B+{xO)) 



< c 



Mc^{B+(xO)) + \Mc{B+{xO)) 



(8.4) 



Proof. We use the a priori interior local Schauder estimate (3.4) for the operator with constant 
coefficients (given by Theorem |3.2[) and the interpolation inequalities for the Holder norms defined 
by the cycloidal metric (Lemma [C.2 ), the method of freezing coefficients as in the proofs of 



[30} Theorem 7.1.1] (elliptic case), [301 Theorem 8.11.1] (parabolic case), and, in particular, [181 
Theorem 3.8] for the parabolic version of our elliptic operator (1.3) to obtain (8.4). □ 



We now generalize Corollary 7.2 to the case of variable coefficients when u has compact support 
in a strip. 

Proposition 8.2 (Higher-order a priori global Schauder estimate for compactly supported func- 
tions on a strip when L has variable coefficients). Let a G (0, 1) and v, Xq, bo, A be positive 
constants and A: S N. Suppose S = M'^""'^ x {0,iy) as in (1.10) and the coefficients a,b,c of L in 
(1.3) belong to Cs'"(S') and obey ( |1.11[ ), ( |1.12 ), and (1.13). Then there are positive constants, 
C = C{k,a,u,d, Xo,bo,A) and 5 = 5(A;, a, i^, d, Aq, ^O) ^) < such that the following holds. If 



u G Cs''^~^'^{S) has compact support in S with diam(suppu) < 6 and u = on diS, then 



and, when c>0 on S, 



< C\\Lu\\^k,C.^gy 



(8.5) 



(8.6) 



Proof Fix x° £ Sn suppu and let L^o denote the operator with constant coefficients a{x^), b{x^), 
c{x^). By applying (7.4) for the operator L^q with constant coefficients, we obtain 



|lt||j-^fc,2 + a^^-j < Co 



\L^gu\\^k.c.f^gy + ||n||c(5) I , 



and hence 



\Lu\\^k,a:^g.^ + 



(L 



(8.7) 



where Cq has the dependencies stated for the constant C in the estimate (7.4). 

For any x^,x'^ € suppn, the cycloidal distance-function bound (2.2) and our hypothesis on 
suppu imply that s{x^,x'^) < \x^ — x^S^I"^ < 6^^'^, for some 6 G (0, z^/2) to be selected later. We 
first consider the case suppn C B^^^y^) for some y^ G OqS. We further restrict to the case A; = 
initially. Observe that 



{L-L 



-a;dtr((a - a(x°))L>\) - {b - 6(x°)) -Du^ic- c(x°))m. 



We consider in turn each of the three terms appearing in our expression for (L — Lxq)u. From 
Definition [2:21 



II (6 - 6(^0)) • DuWc^^s) = Ub- b(.x')) ■ Du\\c^s) + Kb " b{x')) ■ Du]cc 



^It is enough to require It £ C's^°' {B_rg{^'^)) since the estimate trivially holds if ||Lit|||^„j^+ j^oj) and ||it||^(^+ ^^oj) 
are not finite. 

'^This method is also employed in the proof of |26l Theorem 6.2], but Gilbarg and Trudinger employ a family of 
"global" interior Holder norms (which we do not develop in this article) which allows a rearrangement argument. 



REGULARITY OF SOLUTIONS TO DEGENERATE-ELLIPTIC EQUATIONS 



37 



The coefficient bounds (8.14) ensure tliat 



11(6 - 6(:eO)) • DuWc^s) < {\\b\\c{S) + < 2A\\Du\\c^s), 

while the interpolation inequahty ( C.3| ) yields, for some m = ni{d,a) and Ci = Ci{d,a,6) 
(because diam(suppu) = 6) and any e G (0, 1), 



\Du\ 



C{S) 



+ Cie- 



lc{5)> 



u 



C(S)- 



and thus, combining ( |8.8| ) with the preceding inequahty, yields 

11(6 - b{x^)) ■ Du\\q(^§^ < 2eA||M||^2+c<(_g) + CiAe" 

Writing, for x^, x'^ £ S Ci suppii, 

(6(x^) - 6(x°)) • Dujx^) - (6(x^) - 6(x°)) • Dujx^) 



(8.8) 
(8.9) 



we obtain 



[{b-b{x''))-Du]co.^s) < [b]cns) (\\Du\\cis) + 4x\x''nDu]c^^s) 



Since diam(suppn) = 6 and x^,x^ G suppu, by combining the preceding inequality with the 
coefficient bounds (8.14) and the interpolation inequality (8.8), we see that 



[(6 - 6(2;°)) • Du\c^(^s) ^ ^ (^ll^llcJ+°(5) + ""11^110(5) + '5"''^||n||^2+« 



^.10) 



Therefore, by combining (8.9) and (8.10), we obtain 

11(6 - 6(x°)) • DuWc^i^s) < A(3e + 5"/2)||^z||^2+.(g) + 2C7iAe— ||n|b(5)- (8-11) 

An identical analysis, just replacing the coefficient vector 6 by the matrix a, and Du by XdD'^u, 
and the interpolation inequality (C.4) by (C.5), yields 



tr(xd(a - a{x^))D^u)\\co.^s) < + ^"^^)\H\c^+'^^s) + 2CiAe 



lc(S)- 



^.12) 



Similarly, replacing the coefficient vector 6 by the function c, and Du by u, and the interpolation 
inequahty ( [U^ by ( [02] ), yields 

||(c - c(x°))n||c.(s) < A(3e + 



2CiAe' 



We combine (|8.11|), (|8.12D, and (Islsl) to give 



(L-L 



< 3A(3e + (5"/2) 11^11^2+0^^^ + GCiAe- 



\C{S)- 



\c(s)- 



.13) 



We now choose e > such that 9CoAe = 1/4 and choose 6 G (0,i^/2) (which we fix for the 
remainder of the proof) such that 3CoA(5"/^ < 1/4 and combine the preceding inequality with 
( [8^ to give 

for some constant C2 with at most the dependencies stated for C stated in our hypotheses. 
Rearrangement and the maximum principle estimate (Corollary A. 2) for ||ii||(7(5) now give the 
conclusions (8.5) and (8.6) when c > on S" in the case k = 0. 
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Next, suppose k > 1 and let /3 G N'^ be a multi-index with |/3| < k. Because 

D^{vw)= Yl dP'^D^"w- 

l3'+l3"=/3 

for any v,w ^ Cs'"(S'), we may apply the preceding analysis virtually unchanged with v = 
a — a{x^), b — b{x'^), or c — c{x^) and w = XdD'^u, Du, or u, respectively, for each P,(3',I3" G N*^ 
with < k and fi' + f5" = /?. This completes the proof when suppu C B^{y^) for some G d^S. 

Because suppu C Bg/2{x*), for some x* G S, the case dist(j;*, 9o'S') < 5 is covered by our 
analysis for half-balls, with G doS. If dist(x*, ^oS") > S/2, then the operator L is 

strictly elliptic since Xd > 5/2 and |26| Theorem 6.6 and Problem 6.2] imply that 



k,2+a 



< Co \\Lu 



IC?'"(S)) 



+ u 



C{S) 



which is just (8.5). Combining the preceding inequality with the maximum principle estimate 
(Corollary A. 2) for ||tt||c(5) again gives the conclusion (8.6) when c > on S. □ 



Finally, we use Proposition |8.2| to generalize Theorem |7.1| to the case of variable coefficients to 
obtain the following analogue of [H Theorem 1.1.3] (for a related degenerate-parabolic operator 
with constant coefficients) and [26', Corollary 6.3 & Problem 6.1]. 

Theorem 8.3 (Higher-order a priori interior local Schauder estimate when L has variable co- 
efficients). Let a G (0, 1) and let ro, Aq, bo, A be positive constants and let k £ N. Suppose the 
coefficients a,b,c of L in (1.3) belong to Cs'°'{B_^ (x^)), where x^ G dM, and obey (8.2), (8.3), 



and 



+ l|c|Ufc,ar,p+ f n.s < A. 



^.14) 



Then, for any r G (0,ro), there is a positive constant, C = C{k, a,r,ro,d, Xo,bo, A) such that the 
following holds. If u £ Cs'^^°' {Bj^ (x^)) then 



|Lu|Lfc,Qr,B+ ,..n^^ + \\u\\c(B+^{xO)) 



(8.15) 



Proof. We apply an induction argument. When k = 0, the estimate (8.15) follows from Theorem 



8.1 and so we may assume without loss of generality that k > 1. By induction, we may assume 



that the estimate (8.15) holds with constant C = C{1, *) when k is replaced by / G N in the range 
0</<A;-l. 

Let ri := {r + ro)/2 and choose a cutoff function ip G C^{M) such that < 99 < 1 on 



and If = 1 on Bj. 



while supp(/3 C B:^(x^) and note that, for u G Cs''^~^°'{B_^ (x'^)) and thus 



uo :-- 



(fu G Cs''^^°'{S), we have 



C?'2+"(5+(x0)) 



< 



c: 



k,2 + a 



where S 
obtain 



X (0, ro) is the strip as in (1.10), and uq = on diS. By Proposition 

||Mo||(^fc.2+Q(-g^ < Co (^\\Luo\\^k,o,^g^ + ||uo||c(5)) ; 



8.2 



we 



where we use Co to denote the constant C in (8.5), and hence, combining the preceding two 
inequalities, 

ll^llcF+-{B+(xO)) ^ Co (ll^Mllc^."(S+(xO)) + \Mc{B+{xO))) ■ (8-16) 
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Notice that L[ipu) = ipLu + [L, ip\u and that [L, ip] is a first-order partial differential operator, 
[L, ip\u = L{ipu) — ipLu 

= — tT{xdaD'^ ifu)) — b ■ D{ipu) + ip iv{xdaD'^u) + (fb ■ Du, 

and so 

[L, ip]u = - tr{xda{{D'^ip)u + Dip x Du)) - (6 • Dip)u, (8.17) 
where Dip x Du denotes the d x d matrix with entries pxiUxy Observe that 

[L,p]u\\^,,^^^+^^^,^^ + \\pLu 



C^-"(B+(xO))- 



Because of the structure (8.17) of (with factor Xd in the coefficients of t 

derivatives) and the fact that C7,^'°(S+(x°)) C "^'^+"(5+ (x^)) (by Definitions 
we obtain 



le first-order 
and 



2.2 



2.3), 



||[L,<^]n|| 



C^'"(B+(xO)) 



< C||m| 



C: 



k-\,2 + a 



(B+(xO))' 



where C has at most the dependencies stated for the constant in the estimate (8.15). By our 
induction hypothesis, we can apply the local Schauder estimate (8.15) with k replaced by Z = A: 
to give 



1 



C: 



fe-l,2 + c« 



(B+(a:0)) 



< CiWLu 



cr^'"(s+(xO)) 



+ \\u\ 



CiB+ixO)) 



Combining the preceding three bounds with (|8.16|) yields the inequality 

< C 



Ci'^+°'{B+(xO)) 



\Lu\ 



and this is (8.15). 



(8.18) 

□ 



We can now prove the generalization of Corollary 7.2 to the case of variable coefficients. 



Proof of Corollary 1.3. The proof is virtually identical to that of Corollary 7.2 except that we 
replace appeals to Theorems 3.2 and 7.1 (for the case of constant coefficients) by appeals to 
Theorems 8.1 and 8.3 (for the case of variable coefficients). □ 



We can now give the 



X 



Br,{x^)r\\ 



c 



Proof of Theorem \l.l\ Since we can apply Theorem 8.3 to half-balls, -B+i 
when xq E dff, and the standard a priori interior Schauder estimate for strictly elliptic operators 
pel Corollary 6.3 & Problem 6.1] to balls, B 



(£ when xq € 



argument is very similar to the proof of Corollary |7.2[ when the domain, 
8.2. Regularity. We begin with the proof of our global existence result on strips 
Proof of Theorem \1.6{ The proof follows by the method of continuity. We denote 



, the remainder of the 
, is an infinite strip. □ 



-Xd 



Then Corollary 
We consider Lt : 



B.4 



i=l 



dx 



implies that, for any / G Cs'"(5), there is a unique solution u G Cs''^~^"{S). 
1 — t)^Q + tL, for all t G [0, 1]. Given Corollary 1.3 and the existence and 
uniqueness of solutions in '^"''"(5) for the operator .ifg, the method of continuity [26l Theorem 
5.2] applies and gives the result. □ 

We have the following analogue of |26^ Theorem 6.17], albeit with a quite different proof. 
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Theorem 8.4 (Higher-order interio r loc al regularity of solutions when L has variable coefficients). 
Assume the hypotheses of Theorem 8.3 for the operator L in (1.3). If u ^ C'^{B^ (x^)) obeys 



u, Du, XdD'^ueC{B+{x^)) and Lu e C^'"{B+(x")), 



„0^ 



XdD'^u = ondoB^ix^), 



(8.19) 
(8.20) 



then u £ C[ 



k,2+a 



Proof. Let r G (0,ro), and ri := (r + ro)/2, and r2 := (ri + r)/2, so r < ri < r2 < tq. 
Let if G C^{M.) be a cutoff function such that < (/? < 1 on H with (/? = 1 on ^+(x'') and 
supp99 C S+(x°). Denoting S = R'^"^ x (0,ro) as in (1.10) and uq := uip on S, we see that 
uq E 0^(5) is a solution to (1.16), (1.17) with / replaced by 

/o := ipLu + [L, Lp]u on S. 

By hypothesis, (pLu G Cs'"(i?^(x'')), while Lemma 
[L,^]ueCf{S), so 

h G (5), 

while the conditions (8.19) and (8.20) on u imply that uq obeys 

uo, Duq^ x^D'^uq G C{S_) and x^D'^uo = on doS. 



C.3 



and (8.19) and (8.20) ensure that 



Corollary B.4 implies that there is a unique solution V G Cf+"(S') to (fTTel), (fTTTl) with / replaced 



by /o and the maximum principle. Lemma 



A.l 



implies that no = v. Thus, u G C2+"(S+(x°)). 
When A; > 1, we argue by induction and suppose that u G Ct^'^+'^iBlix^)) as our induction 



hypothesis. Bu t the n [L,ip]u G Cs'°'{B^{x^)) by the proof of Theorem 8.3 and so /o G Cs'"(5). 

(S) by the preceding argument for A; = 0, and thus 

□ 



Now Corollary B.4 implies that f G C, 



k,2+a 



u & Cs' °'{B^{x )) since uq = v on S. 
We can now complete the 



Proof of Theorem \1.8[ This is an immediate consequence of Theorem 8.4 and [26', Theorem 6.17] 
since we can apply those regularity results to any half-ball, B:^ (x^) = BrQ{x^) n IH (s ^ when 



' G (9BI, or ball, BrQ{x^) (s when x^ G H, respectively. 
Finally, we complete the proofs of Theorem |1.11| and Corollary |1.13[ 



□ 



Proof of Theorem \l . 1 1\ The argument is very similar to the proof of Corollary B.4, so we just 
highlight the differences. Because / G Cs'"{ff) H Cb{^), we can apply the regularizing procedure 

Tl 

/inC.''"(c7)na 



described in |6l §L11] and [H Theorem 1.11.3] to construct a sequence of functions {/njneN C 
C^(H) such that ^ / in Cs'"(C7) n Ch(^) as n ^ oo, for aU U ff, and 

\\fn\\^^,c.,rr,.<C'\\f\\^k,c„rr., Vn G N, 

where U' <£ U_ and U <£ ^ and C" may depend on U and ?7', and 

||/n|lc(^) < C7||/|b(^), VnGN, 

for some positive constant, C = C{d). 

Let {ttnln gN C C"^(^ ) n C(^ U 5i^) n Cfc(^) be the corresponding (unique) sequence of 
solutions to ( |1.22 ), (1.23), with / replaced by /„, provided by [191 Theorem 1.11]. The maximum 
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principle estimate (Corollary A. 2 for the case cq = and [16, Proposition 2.19 & Theorem 5.3] 
for the case cq > 0) implies that 

for a constant Co depending on the coefficients of A and v when height (^) = v and cq = or 
Co = 1/co when co > and height (^) = oo. The remainder of the argument is now the same as 
the proof of Corollary |B.4[ □ 



Proof of Corollary \ 1. 13\ The conclusion follows immediately from Theorem 1.11 and |19l Corol- 
lary 1.13], since the latter result ensures that u € C{G'). □ 



Appendix A. Maximum principle for degenerate-elliptic operators on domains of 

FINITE HEIGHT 

In this appendix, we prove a comparison principle for operators which include those of the 
form L in (1.3) with c > when the domain, G, is unbounded. Notice that when c does not 
have a uniform positive lower bound, the weak maximum principle Theorem 5.3] does not 
immediately apply when G is unbounded. 

Lemma A.l (Comparison principle on a strip). Let G be a domain of finite height. Ler^' 



Lv := — tr(a-D^t;) — b ■ Dv + cv on 



require that its coefficients, a : 



^ R'^'"^, andb: G- 
a{x) = on 
«, > on G, G R'^.^ 

inf 6^^ > on G, 

e 

tr(a(x)) + (x,6(a;)) < i^(l + Ixp), Vx G 

sup af''^ < oo on 
e 

c > on 

for some •positive constant K. Suppose that u G C^(^) n Cic 
Du, ti{aD^u) G C{G), and 

tT{aD^u) = on d^G . 
If Lu < on G and u < on diG, then u < on G. 

Proof. Define constants bo > and A > by 

A := sup a'^'^ and bo := inf b'^. 

Let o" be a positive constant, to be fixed shortly, and define v G C^ 
by the transformation 

u{x', Xd) = e~'^'''^v{x', Xd), yix', Xd) G 



and c : ^ — )• M obey 



and sup^n < oo, and 



noting that sup^ f < oo since sup^ f < oo and heig^ 
calculation, we find that 



(A.l) 

nCioc(i^) with svipffV < oo 
(A.2) 

< oo by hypothesis. By direct 



Lu 



2aa''^ ]vx, + {c + ab' 



a a 



10 



Note the more general definition of tlie coefficient a{x) in Lemma A.l and Corollary 



A.2 
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We now define coefficients d, b, c of an operator L by 



Lv :-- 



VVx, + CV 



2aa''^]vx. + ic + ab'^-a^a'^'^]v, 



and we notice, by our liypotheses on u and definition of v, tliat Lv < on ^ and t> < on di^. 
Since a = on cJq^, we have 

b"^ = b'^>bo>0 on doff. 
We now choose a := bo/{2A), so that, using c > and a'^'^ < A, 

c = c + ab'^- a'^a'^'^ > abo - cr^A = abo/2 > on ^. 

Then, |! 16| Theorem 5.3] applies to v, and now the conclusion foUows immediately for u also. □ 

Corollary A. 2 (Maximum principle estimate). Let u > and let & ^ x (0, v) be a domain, 



let L be as in Lemma A.l . and let f £ C t,{ff), and g G Cf,{diff). If u obeys the regularity properties 
on G in the hypotheses of Lemma A.l and 



Lu = f on 0', 
u = g on di^, 



then there is a positive constant, C = C{u,bQ, h) with bo, A as in (A.l), such that 

\c{^) + IbllccO]^)) • 



Proof. We define v = e"^''-u as in (A. 2), where a is chosen as in the proof of Lemma A.l Then, 

Lv = f on 1^, 
V = g on di^, 

where / := e'^^'^f on and g := e'^^'^g on dii^. Because c > 6g/(4A) > on ^ from the proof of 
Lemma A.l, we can apply [,16] Proposition 2.19] to give 



1 

< 



C(d) + llffllc(9i^ 



The conclusion follows since Xd G [0, v] for all x & ff and we can take C := e°"^/c = AAe'^'^ /b^. □ 



Appendix B. Existence of solutions for degenerate-elliptic operators with 

constant coefficients on half-spaces and strips 

In this section, we prove existence of smooth solutions to Lu = / on the half-space IHI or on 
'^"^ X (0, ly) as in (1.10), for some u > 0, when the source function / is assumed to be 

"^^^ X [0, z^), respectively, and under the assumption of 



strips S - 

smooth with compact support in HI or in 



Hypothesis 3.1, that the coefficients, a,b,c, of the operator L in (1.3) and so the coefficients, a,b 



of the operator Lq in (1.4) are constant. The method of proof is similar to that of [6, Theorem 
1.1.2] and it is based on taking the Fourier transform in the first {d — l)-variables. The problem 
is then reduced to the study of the Kummer ordinary differential equations whose solutions can 
be expressed in terms of the confluent hypergeometric functions, M and U [TJ §13]. 

We begin by reviewing the properties of the confluent hypergeometric functions which will be 



used in the proofs of Theorems 1.5 and B.3 
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Lemma B.l (Properties of the confluent hypergeometric functions). [T] Let a G C 6e such that 
its real part is positive, K(a) > 0, and let b be a positive constant. Then the following holds, for 
all y > 0. 

(1) Asymptotic behavior as y ^ +oo: 

M{a, b, y) = (l + 0{y-^)) P § 13.1.4], (B.l) 

U{a,b,y) = y-''{l + 0{y-^)) [D §13.1.8]. (B.2) 

(2) Asymptotic behavior as y ^ 0.' 

M(a, 6, 0) = 1 ili§ 13.1.2] and M{a,b,y) = I + 0{y) §13.5.5], (B.3) 

and 

C7(a,6,y) = ^^yi-'' + 0(y^-2) ^^^^3, [1, §13.5.6], 

r(a) 

U{a,b,y) = ^^^^^y^-^ + 0{\\ogy\) ifb = 2, ^1, §13.5.7], 
r(a) 

U{a,b,y) = ^^1~^\ ^''> + 0{1) ^/l<6<2, [1, §13.5.8], 

_ (B.4) 
U{a,b,y) = -^^{\ogy + ^{a)+2^)+0{y\\ogy\) ifb = l, |Ii §13.5.9], 

U{a,b,y) = r(7|a-\) i/0<6<l, [1, § 13.5.10], 

^("'^'0) = ^ ^/0<6<l, m §13.1.2 & 13.1.3], 

1 (1 + a — oj 

where 4^{a) = T'{a)/T(a) and 7 G M is Euler's constant [1, §6.1.3]. 

(3) Differential properties: 

M'{a,b,y) = ^M{a + l,b+l,y) [1, § 13.4.8], (B.5) 
U'{a,b,y) = -aU{a + l,b + l,y) [1, § 13.4.21]. (B.6) 

(4) Recurrence relations: 

(6 _ i)M(a -l,b-l,y) = {b-l- y)M{a, b, y) + yM'{a, b,y) ^ § 13.4.14], (B.7) 
U{a-l,b-l,y) = {l-b + y)U{a,b,y)-yU'{a,b,y) P §13.4.27]. (B.8) 

We can now give the 

Proof of Theorem \1.5\ Uniqueness of the solution, u G C°°{M), follows from the maximum prin- 
ciple [161 Theorem 5.3]. By simple changes of variables described in the proof of [18, Proposition 
A.l], which leave invariant any strip of the form M"^^^ x (0, u), for > 0, we may assume without 
loss of generality that a*-' =6^^ in (1.3) and so the differential operator L has the form 

Lv = -XdVxiXi - b'vxi + ct", G C^(E[), 
where b'^ and c are again positive constants. 
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We adapt the method of proof of [6, Theorem 1.1.2]. We fix / e C^(a). If u G C°°{m) is a 
solution to Lu = / on H, then we expect that its Fourier transform in the x' = {xi, . . . ,Xd-i)- 
variables, 

is a solution, for each ^ G M'^"^, to the ordinary differential equation, 

^rf) + (^c + i y^ik + \i?x}^ u{i- xd) = /(e; xd), (b.o) 

for all Xd G (0,oo), where f{^',Xd) is the Fourier transform of f{x\xd) with respect to x' G M'^. 
We show that the ordinary differential equation (B.9) has a smooth enough solution, tt, in a sense 
to be specified, such that its inverse Fourier transform, 

u{x\xd):=—\^ I u{C;xd)e''^'^dC, (B.IO) 
(27r)'='/^ J^d-i 

is a C°°(]HI) solution to the equation Lu = / on H. 

Defining the function v{^; y), for y = 2\^\xd and each ^ G M*^"^ \ {0}, by 

n(e;xrf) =:e-l«l^''^i-(e;2|C|xrf), G M^-^ \ {0}, Vx^ G M+, (B.ll) 

we see that is a solution to the Kummer ordinary differential equation, 

-yVyy{^;y)-{b-y)vyiC;y) + ai0vi^;y)=9i^;y), VyGM+, (B.12) 

where we denote 

:= — , 



y 



Because > and c > by hypothesis, we see that b > and 3?(a(^)) > when ^ / 0. Since / 
has compact support in H, the function g{(^; •) also has compact support in M+. 

It suffices to study the solutions, v{^; •), to the Kummer equations for ^ G W^^^ \ {0}, and so 
without loss of generality, we will assume in the sequel that ^ 7^ 0. The remainder of the proof 
of Theorem |1.5| is completed in two steps. 



Step 1 (Solution to the Kummer ordinary differential equation). The general solution to the 
Kummer ordinary differential equation ( B.12[ ) can be written in the form v = +vP, where 

v^'i^; y) := ciM{a{0,b; y) + C2U{a{i), b; y), 

v^{i-y) ■.= -M{a{0,b;y) H g{i; z)^f%^dz - U {a{i)My) f 9{i;z)^P§^^dz, 

Jy W{a{^),b;z) Jq W{a{^),b;z) 

with ci,C2 G M, and 

WiaiO,b;y) :=- ^ffi''f , G R'^-^ {0}, Vy G M+, (B.14) 

is the Wronskian of the Kummer function, M{a{S,),b,y), and the Tricomi function, U{a{(,),b,y), 
P § 13.1.22]. We want to find a solution, v G C°^(t+), to (IbIIJ). 
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From (B.l ), we see that the function M{a{^), b; y) is unbounded as y tends to +00, and so we 
choose the constant ci = 0, because we only consider bounded solutions. At y = 0, we obtain 
from (B.6) and (B.4) that U'{a{S), b, y) is unbounded, since 6 > 0, and so we choose the constant 
C2 = 0, because we only consider solutions to the Kummer equation which are smooth on M+. 
Thus, we obtain 



v{i;y) 



-M{a{0,b-y) 
-U{a{0,b;y) 



' ^ U{a{0,b;z) ^ 
W{a[Q,b]z) 



^^^'"V(a(0,6;z)^"- 



y 



(B.15) 



Given v defined as above, and u defined as in (B.ll ), we will prove the following properties of the 
solution, u, to verify that u defined by (B.IO) is a C°°(]HI) solution to Lu = / on H, as asserted 
by Theorem |1.5[ 



Lemma B.2 (Properties of u). If f ^ C[ 

following properties. 

(1) For all i G M^-^ \ {0}, we have 



, then the function u defined by (B.ll) has the 



lim u{i;xd) 



0. 



(B.16) 



(2) The function u{^; •) belongs to 

(3) The function n(^; •) obeys 

1 



for all ^ G 



\{0}. 



C y>0 



\ {0}, Vxrf G 



(B.17) 



where c is the zeroth-order coefficient of L in (1.3). 

(4) The function u{-;xd) decays faster than any polynomial in ^, for all Xd G ^ 

(5) The functions u decay faster than any polynomial in ^, for all G N. 



Step 2 (Existence of a solution, u G C°°{M), to Lu = / on H). From Lemma B.2 Items ([2]) and 



([4]), we see that the function u defined by ( |B.10 ) has an arbitr ary number of derivatives in the 
first (c? — l)-variables which are continuous on H. From Lemma 



B.2 



Item ([2]) , we see that u also 
admits an arbitrary number of derivatives in the x^-variable, and they are continuous on H. Now 
we consider D^'^^'^u, for j3d G N, which satisfies 



D^^''^u{x\ Xd) := / D^^'^u{i- xrf)e*^'« d^, Vx^ G 



By Lemma 



B.2 



Item ([5]), the function D^'^^'^u decays faster than any polynomial in ^, and so 
D^D^dedy^ exists and is continuous on H, for all A G N'^ with = 0. Thus, u belongs to C°°(BI). 
Since u solves (B.9), we find that u solves Lu = / on IH by taking the inverse Fourier transform 
of u{C;xd) in ^ G M'^"^ From {W, Theorem 5.3], it follows that u is the unique C°°(]HI) solution 
to Lu = / on H. 



Aside from the proof of Lemma B.2 given below, this completes the proof of Theorem 1.5 □ 



It remains to prove Lemma B.2 



Proof of Lemma B.2. We organize the proof into several steps. 



46 



P. FEEHAN AND C. POP 



Step 1 (Proof of Item ([Ij)). First, we verify that the function v defined in Equation (B.15) is 
well-defined. We write v = + , where we set 



y 



Mia{0,b; z) 



dz 

Recall that •) has compact support in M+, and so for the function we only need to 

verify that it is continuous up to y = 0. From the property (B.3) in Lemma |B.1| we know that 
M(o(^), 6; y) is continuous in ?/ up to y = with M(a(^), 6; 0) = 1. Identities (B.4) and definition 
(B.14) of the Wronskian imply that 

U{a{i)My) 
W{a{0,b;y) 



'max{y, y2(b-2)}^ if 6 > 2, 

ylogy, if 6 = 2, 

y, if 1 < 6 < 2, 

ly\ ifO<6<l, 



and so this function is integrable near y = 0. Since g{S,; •) has compact support in 
that vHC, •) G C{R+) and 

lim v^{C;y) = 0, G M'^"^ \ {0}. 

y/'oo 

Next, we consider the behavior of the function •). Near y 
definition (B.14) of the Wronskian yield 



we see 



0, the property (B.3) and 



M(a(0,6;y) 



y • 



0, we find that the limit of 



W{a{C),b;y) 

Combining this result the asymptotic behavior ( B.4[ ) of U as y 
v'^iC'i y), y tends to 0, exists. The limit of the integral, 

W{a{Q,b;z) 

as y ^ CO obviously exists because the function g{£,;-) has compact support in ]R_|_. Moreover, 
using the asymptotic behavior (B.2) of U{a{^),b;y) as y — )• +oo, we obtain 

lim v'^(C;y) = 0, G M'^"^ \ {0}. 

y/'oo 



Since v = + v'^, we obtain the limit property (B.16) for n as y — t- +oo using (B.ll). 

Step 2 (Proof of Item ([2])). The argument employed in Step[l]shows that u{(,; •) G C(M+), for all 
C G M"'"^ \ {0}. Next, we want to show that u{^; •) e C(M+), for ah A; G N and ^ G R'^-^ \ {0}, 



but for this it is suffices to show that D'^v{C; •) G C(M+), for all A; G N, by l\B.ll\) . 
We first consider the case k = 1. A direct calculation shows that 



Vy{C;y) = -My{a{0,b;y) 



g[C;z) „^ — rdz 



Uyia{0,b;y) / g{^;z 



W{a{0,b;z) 
M(a{0,b;z) 
WiaiO,b;z) 



dz. 
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Using identities (B.5) and (B.6), we obtain 



<0 



POO 

M(a(e) + l,6+l;y) / g{C; z) 
Jy 



WiaiO,b;z) 



dz 



+ amiaiO + l,b+l;y) r g{^- z) ^j^ig'?'"! dz, 



(B.18) 







W{a{i),h-z) 



and the same argument as in the beginning of the proof of Lemma [B.2[ gives us immediately that 
Vy{i] •) e C7(M+). Hence, v{i; •) G C^(R +), for all ^ G \ {0}. 

We next show that Vy{£^; •) in (B.18) is the unique C^(M+) solution to the Kummer equation, 

-ywyy{i]y) - (6 + 1 - y)wy{i] y) + {a{i) + l)w{i; y) = gyi^; y), Vy G M+. 

Our goal is to show that Vy = w, where we define 



/•oo 

w{^- y) := -M(a(e) + 1, 6 + 1; y) / z 



U{a{i) + l,h + I- z) 

V(a(o + i,6+i;z; 



dz 



, , . , X .M{a{C) + l,h+l:z) 



for y € M+, ^ G M \ {0}. Integrating by parts in the expression of we obtain 
w{i-y) = M{a{i) + l,b+l-y) j g{^-z)^^^ ^a{i) + l,h + I- z) dz 

+ C/(a(e) + l,6+l;y) g{e, z) '-{a{0 + hb + 1; z) dz. 



The expression of Vy in (B.18) coincides with that of w if 

(a(6 + 1,^ + 1;^) = 

-(a(e) + 1, & + 1; z) = -a{i) — {a{i)M z). 



w 



But the preceding two identities follow from the definition of the Wronskian, W , in (B.14), from 
(B.5) and (B.6), and from the recursion relations (B.7) and (B.8). Hence, the function Vy is the 
unique C^(M+) solution to the corresponding Kummer equation, which obviously implies that 

Vyy £ C(M+). 

Inductively, it follows that, for any A; G N, the derivative DyV exists and is the unique C^(R+) 
solution to the Kummer equation, 

-yiD'yV)yyiC; y)-{b + k- y){Dlv)y{i- y) + (a(e) + k)Dlv{i; y) = D';g{C; y), Vy G M+. 



Thus, viC; •) G C"^(M+), and so u{£,; •) G C"^(M+) by ( [bTT| ), and D''u{C; ■) satisfies the ordinary 
differential equation, for all Xd G (0, oo) and A; G N, 

- XdiD^^,u)^,^,i^; Xd)-{b + k- 2\C\xd)iD^^^u)^,i^; Xd) + (a(0 + k)D^X^; Xd) 



(B.19) 

Notice that the right-hand side in the preceding equation is a function with compact support in 
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Step 3 (Proof of Items ^ and We adapt the argument of [61 Theorem 1.1.2]. We fix ^ / 0. 
We write u{^; Xd) = p{i\ Xd) + iq{S,; Xd) and /(^; Xd) = g(^; Xd) + ih{i; Xd). Then, equation ( [R9 ) 
becomes 



-XdPxaxai^; Xd) - h'^PxM'^ Xd) + (c + Xd\i\'^)p{i] Xd) - hiq{i; Xd) = g{i; Xd), 
-XdQxaxaiCxd) - h'^qxM'^Xd) + {c + Xd\i?)q{i;xd) + b^p{C;xd) = h{^;Xd), 



where 6^ denotes the inner product of (6^, . . . , ^) with ^. Defining F{£^; Xd) 
p'^iC'i Xd) + q'^iCi Xd), we obtain (where now we omit the (^; Xrf)-variables) 

= '^ppxd + '^qqxa, 

Px^xd = ^pI. + 2pp^^^^ + 2ql + 2qq^^^^, 



\u{C;Xd)\'' 



which gives us 

XdFx.x, + b^F,, - 2cF 



= 2p{xdPxaxa + b'^Pxa - cp) + 2q{xdqxaxa + b'^qxa - cq) + 2xd{pl^ + (1: 

> 2p {-g + XdlePj? - b^q) + 2q [-h + Xd\S,\^q + b^p) 

> -2pg - 2qh 

> -cF + 

c \ 

where we recall that c > by hypothesis, and so it follows that 

XdFxaxai^'^Xd) + b'^Fa:a{^;xd) -cF{(;xd) >-- sup \f{^;xd) 



^Xd' 



\/xd G 



Now let 



Then, 



G{i-Xd) ■.= F{i-Xd) 



sup \f{i;xd) 



Xd& 



XdGx^XdiC; Xd) + b'^G^^iC; Xd) - cG{i; Xd) > 0, 
lima;^/ioo G{^; Xd) < 0, 



where we used (B.16) to determine the behavior of G{S^;xd) as Xd approaches c«. Therefore, the 
function G(^; x^) is bounded, and the maximum principle [16, Theorem 5.3] then implies that 
G{C, Xd) < 0, for all Xd G M+, and so 



u{^;xd)\^ < \ sup |/(^;y)P, Vx^ G 



^+1 



which is equivalent to (B.17). 

Since / belongs to Cq°(]HI), the function sup |/(^; y)| decays faster than any polynomial in 



^ by |23| Theorem 8.22 (e)]. Therefore, from (B.17) we see that the function u also decays faster 
than any polynomial in ^. 



Step 4 (Proof of Item ([5])). By (B.19) and the fact that the right-hand side in (B.19) is a function 
with compact support in M+, we see that the preceding steps can be applied toD^^u instead of 
u, for all fc G N. Therefore, we obtain that the functions F>^^u decay faster than any polynomial 
in ^, for all A; G N. 
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This completes the proof of Lemma B.2 □ 
We now prove the existence and uniqueness of smooth solutions on strips in the half-space. 



diS 



-1 



X (0, z^), so that doS 



pd-l 



X {0} and 



We fix 1/ > and recall from (1.10) that S = 

'^^^ X {z^}. We have the following elliptic analogue of [6l Theorem 1.1.2] in the parabolic 



case, but for finite-height strips rather than the half-space. 

Theorem B.3 (Existence and uniqueness of a C°°{S) solution on a strip when L has constant 
coefficients). Let L be an operator of the form (1.3) and require that the coefficients, a,b,c, are 
constant with b'^ > and c > 0. Then, for any function, f e C^{S), there is a unique solution, 
u £ C°^{S), to 

(Lu = f on S, 
u = on diS. 



(B.20) 



Proof. The method of proof is the same as that of Theorem 1.5, so we only highlight the main 
diffe rences. Uniqueness of the solution, u G C°°{S), follows from the maximum principle, Lemma 
a1| for L. By analogy with ^B^, for each ^ G M'^-^ \ {0}, we construct the function u{^; 



to 



be the unique solution in C°°([0, to 



-XdUxaxM'^^d) - b'^u^A^^Xd) + c + i^b^ik + \i?Xd u{i;xd) = f{^;xd), Vx^ E (0,i/), 



k=l 



by defining the new function, v{^,; •), by (B.ll) and proving that v{^; •) is the unique solution in 
C°°([0, xJJ]) to the Kummer equation, 

-yVyyiC; y)-ib- y)vyi^; y) + ai^vi^; y) = g{^; y), Vy e (0, 2\C\u), 
for each ^ G M'^~^ \ {0}, where the coefficients b and a(^), and the function g are defined in the 



same way as in (B.13). The arguments employed in the proof of Theorem B.3 show now that the 



unique solution in C°°{S) to the preceding ordinary differential equation is given by 
v{C; y) := CM{a{0,b, y) - M{a{i)M y) 



g{i;z) „^, . — -dz 



y 



W{a{i)Mz) 



where the constant C is chosen such that the boundary condition, ^(,^,21^12^) = 0, is satisfied. 
The only remaining modification that we need lies in Step [3] of the proof of Lemma B.2 The 
reason why this part of the proof does not adapt immediately is because we used the fact that 



the zeroth-order coefficient, c, of L in (1.3) is strictly positive to derive (B.17), while now we 
assume c > 0. To circumvent this issue, we apply the method of proof of Step |3] of Lemma B.2 
not to F, but to e~'^^''F, where we choose the positive constant, a, small enough. Notice that 
this is the same as the approach we employed in the proof of Lemma |A.1| to overcome the fact 
that c = 0. □ 



Corollary B.4 (Existence and uniqueness of a Cs solution on a strip when L has constant 



coefficients). Let a G (0,1) and /c G N. Let L be an operator as in (1.3) and require that the 



coefficients, a,b,c, are constant with b'^ > and c > 0. If f £ Cs'"(5), then there is a unique 



solution u G C^''^~^'^{S) to the boundary problem (B.20). 
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A.l 



Proo f. Uniqueness of the solution, u E Cs''^^"(S'), follows from the maximum principle, Lemma 
for L sinc e any u £ Cs'^^°(5') has the property t hat D u and x^D'^u are continuous on S by 
and that XdD'^u = on OqS by Lemma C.l Let {/njneN C Cq^{S) be a sequence 
such that fn^f in Cs'^~''"(5') as n — )• cxd and 



Definition 



2.3 



Such a sequence can be constructed using j6l Theorem L11.3]. Let u„ S C°^{S) be the unique 
solution to (B.20), with / replaced by /„, given by Theorem B.3 In particular, each solution 



satisfies the global Schauder estimate (1.15) which, when combined with the preceding inequality, 
gives 

\\un\\^k,2+c>^g-^ < C||/||^fe,c«^^^, Vn e N. 

By applying the Arzela-Ascoli Theorem, we can extract a subsequence, which we continue to 
denote by {unlneN; which converges in Cs''^~^°'{S) to a limit function u G Cs'^^"(5') as n — )■ 
oo. Since {fn}neN, {un}neN, {Dun}neN and {X(iD'^Un}neN also converge uniformly on compact 
subsets of S to /, u, Du and x^D'^u, respectively, as n — )• oo, we see that u solves (B.20). □ 



Appendix C. Interpolation inequalities and boundary properties of functions in 

WEIGHTED Holder spaces 

A parabolic version of following result is included in [6^ Proposition 1.12.1] when d = 2 and 
proved in [18] when d> 2 for parabolic weighted Holder spaces. For completeness, we restate the 
result here for the elliptic weighted Holder spaces used in this article. 

Lemma C.l (Boundary properties of functions in weighted Holder spaces). [18^ Lemma 3.1] // 

u G Cf+^d) then, for all x^ £ dW, 



lim XfiD u{x) = 0. 



(C.l) 



In [18] , we also proved the following interpolation inequalities parabolic weighted Holder spaces 
analogous to those for standard parabolic Holder spaces [30l|32]. For completeness, we restate 
these interpolation inequalities below for elliptic weighted Holder spaces, analogous to those for 
standard elliptic Holder spaces in [26l Lemmas 6.32 & 6.35], t30» Theorem 3.2.1]. 



Lemma C.2 (Interpolation inequalities for weighted Holder spaces). [181, Lemma 3.2] Let a € 
(0, 1) and rg > 0. Then there are positive constants, m = m{d, a) and C = C{d, rg, a), such that 
the following holds. If u e Cf+°(S+ (x°)), where x^ G dM, and e G (0, 1), then 



+ Ce-"'\\u\ 



\c{B+{xO)) 

We add here the following 



DuWr^r 1:,+ ,^0-\\ — ^\W\\Cs'^°'{B^i^{x'-')) ~^ ""ll^llr'l'R+ I'^rO 



C(B+(xO))^ 

c{B+ixO)y 

C{B+{xO)y 



\u\ 



C(S+(xO)) 



(C.2) 
(C.3) 
(C.4) 
(C.5) 



Lemma C.3 (Holder continuity for XdDu). Let r > 0, and assume that u G C'^(Br) is such that 
Du and XdD^u belong to C{B^). Then, x^Du G Cf{B^). 
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Proof. For this we only need to show that for any x^,x^ G i?,^ such that aU their coordinates 
coincide, except for the i-th one, we have 

\x\Du{x^) - xlDu{x^)\ 

s{x^,x'^)" ~ 

for some positive constant, C. We show this for the case i = d, and all the other cases, i = 
1, . . . ,d — 1, follow by a similar argument. We have 

\x\Du{x^)-xjDu{x^)\ kd-^ll , o lDnjx^) - Du{x^)\ 

s{x\x^y - six\x^)"^ ^ s{x\x^)^ 

< ( \\Du\\r.,fi+ , + x^,\D\(x^y'^ '"^^ " 



where G i?r-3 is a point on the line connecting x^ and x^, and we apply the Mean Value 
Theorem. We may assume without loss of generality that < x^, and because x^ > x^, we have 
that x^|L'^n(x^)| < ||xrfD^n||^^^+ ^. Using the definition (2.1) of the cycloidal metric, we obtain 

\xlDu{x'^) - xlDu{x^)\ f 2 II \ l-a/2 

Therefore, x^Du belongs to Cf{B^^), for all a G (0,1). This completes the proof of Lemma 
EJ □ 
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